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CHAPTER  1 


INTRODUCTION 


1. 1.  Preview 

This  study  considers  optimal  decision  strategies  with  regard  to 
capacity  expansion  in  a stochastic  demand  environment.  Chapter  1 is  an 
introduction  to  the  topic  at  hand  and  provides  the  preliminaries  upon 
which  later  model  construction  is  based.  As  indicated  in  the  next  section, 
there  are  two  types  of  facilities  which  must  be  considered  in  a capacity 
expansion  model:  permanent  and  temporary.  With  regard  to  temporary 

facilities,  a further  classification  is  necessary  according  to  whether 
or  not  the  temporary  facilities  are  modular. 

Model  I,  presented  in  Chapter  2,  treats  the  case  where  temporary 
facilities  are  non-modular.  Model  II,  presented  in  [17],  treats 
the  case  where  temporary  facilities  are  modular.  For  both  cases,  it  is 
shown  that  the  determination  of  optimal  expansion  sizes  can  be  expressed 
as  a single-variable,  near ly- unconstrained,  minimization  problem  of  a 
very  particular  form.  The  solution  to  this  problem  is  treated  in  118], 
where  a revised  model  is  also  introduced  to  demonstrate  a number  of 
generalizations  that  are  possible  in  both  Models  I and  II. 
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.2.  Permanent  and  Temporary  Facilities 


When  one  considers  the  growth  in  capacity  of  a singular  system, 
such  as  a plant,  pipeline,  superhighway  or  school,  a common  phenomenon 
often  arises.  Specifically,  when  demand  first  exceeds  the  initial  capacity 
of  the  system,  a permanent  capacity  expansion  is  not  immediately  under- 
taken. Instead,  temporary  operating  measures  are  instigated  in  order  to 
satisfy  excess  demand  over  the  short-run.  In  some  cases,  the  temporary 
measure  taken  may  be  simply  to  backlog  excess  demand.  In  other  cases, 
where  backlogging  is  undesirable,  the  present  system  may  be  overloaded 
beyond  its  planned  capacity  in  order  to  accommodate  excess  demand.  In 
still  other  cases,  where  backlogging  is  undesirable  and  overloading  is 
infeasible,  temporary  capacity  may  be  rented  from  outside  sources.  In  order 
to  distinguish  between  the  normal  accommodation  of  demand  and  the  emergency 
accommodation  of  excess  demand  through  temporary  measures,  facilities 
will  be  categorized  here  as  either  permanent  or  temporary.  That  is, 
"temporary  facilities"  will  denote  the  measures  utilized  to  accommodate 
excess  demand  over  the  short-run  prior  to  the  implementation  of  a full- 
fledged  expansion  of  permanent  facilities. 

The  use  of  temporary  facilities  usually  involves  some  sort  of  cost 
penalty  over  that  of  permanent  facilities;  otherwise,  one  might  well  ask 
which  facilities  are  really  temporary.  Thus,  the  economic  advantages  of 
temporary  facility  usage  might,  on  the  surface,  seem  questionable.  There 
are,  however,  two  principal  factors  which  promote  the  desirability  of 
temporary  facility  usage.  Firstly,  there  is  the  problem  of  uncertainty : 
when  demand  first  exceeds  the  permanent  capacity,  is  the  excess  a short-term 
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peak  or  a long-term  trend?  If  the  time  interval  of  excess  is  to  be 
of  small  duration,  then  an  expansion  of  permanent  facilities  is 
undesirable,  since  overhead  costs  on  unused  capacity  can  often  be  very 
substantial.  Secondly,  there  is  the  problem  of  capitalization:  permanent 
capacity  expansion  costs  often  exhibit  significant  economies  of  scale, 
which  in  turn  dictate  substantial  expansion  sizes  requiring  large  capital 
expenditures.  Even  if  sustained  demand  growth  is  certain,  bankruptcy  is 
possible  over  the  short-run  if  sufficient  revenue  cannot  be  generated  from 
new  facilities  to  meet  capitalization  costs.  Furthermore,  capitalization 
itself  may  prove  difficult  until  excess  demand  is  of  sufficient  magnitude 
to  convince  investors  that  permanent  capacity  expansion  is  warranted. 

While  each  of  the  above  problems  can  in  itself  promote  the  use  of 
temporary  facilities,  the  two  taken  together  generate  an  even  stronger 
case  for  temporary  facility  usage  under  an  expected  discounted  cost 
criterion.  Figure  1.1  provides  some  illustration.  The  demand  peaks 
at  times  t^  and  are  short-lived;  any  large  permanent  expansion 

at  these  times  will  result  in  a great  deal  of  wasted  capacity  over  a 
significant  portion  of  the  time  interval  shown.  Temporary  measures 
can  be  used  to  accommodate  these  excesses.  Similarly,  at  time  t^,  the 
long-term  prospects  for  demand  are  unclear  and  temporary  facilities  can 
be  used.  However,  by  time  t^,  the  excess  demand  has  grown  to  significant 

size  so  as  to  somehow  "warrant"  a permanent  expansion  (at  which  time,  the 
temporary  facilities  in  use  would  be  discontinued).  If  temporary  facilities 
are  not  used,  then  a permanent  expansion  must  be  undertaken  at  time  t^. 


Thus,  the  use  of  temporary  facilities  forestalls  the  permanent  expansion 


FIGURE  1.1:  Use  Temporary  Facilities  to  Forestall  Permanent  Expansion 


for  (t^-t^)  time  units.  If  costs  are  discounted,  then  the  discount 
savings  alone,  accruing  from  the  prolongment  of  the  permanent  expansion 
capital  outlay,  can  easily  exceed  any  cost  penalties  resulting  from  the 
use  of  temporary  facilities  at  times  t^  t^  and  t,,. 

In  order  to  make  these  ideas  more  precise,  let  k denote  the 
spares  level:  the  difference  between  present  permanent  capacity  and 

demand.  Then  -k  denotes  the  excess  demand  whenever  demand  exceeds 
permanent  capacity.  When  k is  nonnegative,  the  permanent  facilities 
suffice  to  serve  all  demand.  However,  when  k is  negative,  excess 
demand  exists  and  temporary  facilities  must  be  used.  Let  K denote  the 
limit  on  temporary  facility  usage.  Whenever  k = -K  ( that  is,  excess 
demand  equals  K)  and  a new  unit  of  demand  growth  occurs,  a permanent 
expansion  of  size  X+l  (X  > K)  is  required.  When  permanent  expansion 
occurs,  the  temporary  facilities  in  use  are  discontinued  and  the  spares 
level  k increases  to  X-K  > 0.  If  demand  consists  of  deterministic 
constant  growth,  then  the  use  of  this  type  of  recurrent  (X,K)  expansion 
policy  results  in  a saw-tooth  pattern  for  the  spares  level,  as  illustrated 
in  Figure  1.2. 

The  value  K serves  as  a "trigger",  forcing  a permanent  expansion 
whenever  a new  unit  of  demand  growth  occurs  while  k = -K.  Thus,  K+l 
represents  the  point  at  which  excess  demand  warrants  a permanent  expansion. 
In  cases  of  overloading,  the  value  K may  be  determined  based  on  physical 
constraints  governing  the  degree  of  overloading  that  can  safely  be  incurred. 
In  other  cases,  K may  be  determined  according  to  a judgment  concerning 
the  willingness  of  investors  to  participate  in  capitalization  and  the 
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prospects  of  generating  sufficient  revenue  to  meet  capitalization  costs, 
given  that  excess  demand  has  reached  the  value  K.  In  such  cases  where 
K is  predetermined,  the  parameter  of  interest  will  be  X*(X)+1  the 
optimal  permanent  expansion  size  to  be  utilized  in  conjunction  with  the 
temporary  facilities  usage  limit  K.  In  other  cases,  a value  K*  will 
be  determined  based  on  a cost  minimization  over  all  feasible  operating 
policies  (X,K) . Given  a procedure  for  determining  X*(X),  this  last 
problem  reduces  to  that  of  a cost  minimization  over  all  feasible  operating 
policies  (X*(K) , K) . 

To  summarize,  the  development  of  a capacity  expansion  policy 
typically  involves  at  least  the  determination  of  two  parameters:  the 

permanent  expansion  size  X+l  and  the  temporary  facilities  usage  limit 
X.  The  value  X answers  the  question  of  when  to  add  and  the  value  X 
answers  the  question  of  how  much  to  add. 


1.3.  Examples 

The  use  of  temporary  facilities  to  forestall  permanent  expansion 
is  prevalent  in  many  sectors  at  many  different  levels  of  the  economy. 

The  examples  below  illustrate  the  principal  types  of  temporary  facilities 
utilized. 


Example  1.1.  Backlogging.  In  these  cases,  excess  demand  is  simply  held 
in  an  unsatisfied  state  until  such  time  that  permanent  facilities  become 
available,  either  through  permanent  expansion  or  through  a decrease  in 
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the  demand  level  actually  being  served.  A backlog  penalty  is  charged, 
per  unit- time,  for  each  demand  unit  backlogged.  This  penalty  includes 
the  cost  of  maintaining  the  backlog  and  also,  in  the  case  of  continuous 
service  systems  (e.g.,  utility  companies),  the  revenues  lost  during  the 
period  of  backlog. 

Example  1.2.  Overloading.  In  these  instances,  the  existing  permanent 
facilities  are  overloaded  beyond  their  design  capacity.  Overloading 
usually  places  additional  strain  on  the  entire  system  and  this  increased 
burden  must  be  accounted  for  in  the  form  of  cost  penalties.  A specific 
example  of  this  type  is  that  of  a power  plant  which  is  utilized  beyond 
its  maximal  efficiency  point  in  order  to  accommodate  greater  demand.  In 
this  case,  the  resulting  loss  in  efficiency  must  be  translated  into  a 
monetary  overloading  cost  penalty.  Overloading  cannot  be  utilized 
indefinitely  whenever  demand  growth  persists;  eventually,  a point  will 
be  reached  where  additional  demand  cannot  be  safely  served.  Thus,  an 
upper  bound  for  the  parameter  K is  dictated  by  physical  considerations 
whenever  overloading  is  in  < -curs. 

Example  1.3.  Jobletting.  In  this  case,  outside  sources  are  utilized  to 
increase  overall  capacity.  A specific  example  of  this  kind  is  the  case  of 

a corporate  division  which  has  outgrown  its  in-house  computer  facilities. 
Additional  computer  time  may  be  bought  from  either  an  outside  supplier 
specializing  in  computer  services  or  possibly  from  another  division  within 
the  same  parent  enterprise.  This  specific  example  is  an  instance  of 
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non-modular  temporary  facilities,  (as  are  Examples  1.1  and  1.2),  because 
the  division  can  typically  rent  outside  computer  time  precisely  equal  to 
the  demand  excess.  This  is  also  a case  where  permanent  expansion  may 
in  fact  mean  replacement . If  the  existing  division  computer  configuration 
is  replaced  by  more- sophisticated  machinery,  then  the  expansion  cost  is 
the  price  of  the  new  configuration  less  any  salvage  revenue  generated  by 
sale  of  the  old  configuration.  Similarly,  the  expansion  size  is  the 
difference  in  computing  power  between  the  new  and  old  configurations. 


Example  1.4.  Modular  Temporary  Facilities.  This  example  can  be  viewed 
as  a specific  case  of  jobletting  where  temporary  facilities  can  be  rented 
from  outside  sources  in  a fixed  increment  size.  An  instance  of  this  type 
is  the  use  of  mobile  "barracks"  facilities  (e.g.,  house  trailers)  as 
auxiliary  classrooms  to  alleviate  conjestion  in  overcrowded  schools. 

Another  instance  of  this  type  is  the  use  of  "pair-gain"  devices 
to  augment  the  capacity  of  over- subscribed  telephone  cables.  To  illustrate, 
consider  Figure  1.3.  In  order  to  simplify  the  illustration,  consider  only 
the  one-way  communication  from  a subscriber  telephone  to  the  switching 
equipment  at  the  subscriber's  local  telephone  office.  As  shown  in  Figure 
1.3a,  this  communication  is  normally  implemented  by  the  dedication  of  a 
single  wire- strand  from  the  subscriber  telephone  to  the  office.  This 
single  wire  is  typically  one  strand  of  a telephone  cable  between  the 
subscriber's  neighborhood  and  the  office.  Suppose  that  the  demand  for 
telephones  in  this  neighborhood  increases  to  a point  at  which  all  wire- 
strands  of  the  existing  cable  are  used.  The  telephone  company,  being 
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required  to  provide  service  to  all  soliciting  subscribers,  must  somehow 
increase  the  cable  capacity  in  order  to  satisfy  new  demands.  One  obvious 
possibility  is  the  immediate  placement  of  a new  cable  (i.e.,  an  immediate 
permanent  expansion) . Another  possibility  is  the  implementation  of 
multiple-party  service;  however,  this  alternative  is  usually  unexceptable 
to  consumers  and  is,  in  fact,  not  used  in  many  areas.  Another  alternative 
to  the  placement  of  an  additional  cable  is  the  utilization  of  pair-gain 
devices  on  the  existing  cable.  As  illustrated  in  Figure  1.3b,  a pair-gain 
device  consists  of  electronics  that  permit  the  transmission  of  a number 
of  simultaneous  conversations  over  a single  wire- strand  by  encoding  the 
conversations  at  one  end  and  decoding  the  conversations  at  the  other  end. 

As  shown  in  the  figure,  these  devices  are  typically  available  in  a fixed 
modular  size  allowing  the  transmission  of  L additional  conversations 
over  a single  wire- strand;  the  "pair-gain  ratio"  of  each  such  module  is 
said  to  be  "L  to  1". 

A cost  characteristic  which  frequently  distinguishes  modular 
temporary  facilities  from  non-modular  temporary  facilities  is  the  existence 
of  instaneous  charges  over  and  above  the  normal  rental  costs.  Specifically, 
an  installation  charge  is  typically  incurred  whenever  an  additional 
module  is  engaged  and  similarly,  a removal  charge  is  incurred  whenever 
a module  is  returned  to  the  outside  supplier.  Thus,  in  the  case  of 
modular  temporary  facilities,  there  can  be  a sub- optimization  problem 
with  regard  to  how  modules  should  be  engaged  and  returned  in  order  to 
avoid  excessive  installation  and  removal  charges. 


1.4.  The  Poisson  Demand  Model 

In  order  to  account  for  the  promotion  of  temporary  facility  usage 
due  to  uncertainty,  demand  will  be  recognized  as  being  stochastic  in 
nature.  Spec ifically,  it  will  be  assumed  that  demand  increases  ("arrivals") 
are  characterized  by  a Poisson  process  [15]  at  rate  A^  > C.  Similarly,  it 
will  be  assumed  that  demand  decreases  ("departures")  are  characterized  by 
a Poisson  process  at  rate  Ag  > 0.  Furthermore,  the  arrival  and  departure 
processes  will  be  presumed  to  be  independent  of  each  other  and  also 
independent  of  the  expansion  policy  chosen. 

Given  the  above  demand  characterization,  define  an  "event"  as 
either  an  arrival  or  a departure.  Then  an  equivalent,  and  more  convenient, 
characterization  is  as  follows:  events  occur  according  to  a Poisson 

process  at  rate  A = A^  + Agj  the  probability  that  an  event  is  an  arrival 
equals  p = A^/Aj  and  the  probability  that  an  event  is  a departure  equals 
q = Ag/A  (q  = 1-p)  [15].  That  is,  the  demand  process  treated  here  can 

be  viewed  as  a random  walk  with  exponential  inter- event  times. 

It  should  be  noted  that  the  above  demand  characterization  is 
pointed  toward  systems  providing  a homogeneous  continuous  service  to  a 
fairly  stable  clientele,  such  as  a power  plant  or  a telephone  service  area. 

It  is  really  not  intended  to  model  actual  clientele  movement  for  cases 
where  customers  are  rapidly  flowing  in  and  out,  such  as  a job  shop.  In 
these  cases,  it  is  suggested  that  "customer  demand"  be  viewed  as  the 
maximum  usage  level  required  of  facilities  during  short  intervals  of 
time  (e.g.,  weeks  or  months). 
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It  should  also  be  noted  that  the  demand  process  is  assumed  time- 
stationary. For  this  reason,  the  system  being  modelled  must  be  in  the 
midst  of  a relatively  stable  period  of  either  sustained  growth  or  sustained 
negative  growth. 

Finally,  notice  that  the  arrival  and  departure  rates  are  presumed 
constant  with  regard  to  the  demand  level.  From  a practical  point  of  view, 
this  is  convenient  since  only  two  estimates  --  a single  arrival  rate  and 
a single  departure  rate  --  are  necessary  to  implement  the  model.  However, 
in  most  realistic  situations,  it  appears  that  the  arrival  and  departure 
rates  will  vary  according  to  the  system  demand  level.  In  recognition  of 
this  fact,  the  revised  model  of  [18]  permits  these  rates  to  vary. 


1.5 . Costs 

The  optimization  criterion  used  here  will  be  the  minimization  of 
all  future  expected  discounted  costs;  r will  denote  the  continuous  dis- 
count rate,  0 < r < 1.  Three  types  of  costs  will  be  allowed:  permanent 

expansion  costs,  permanent  facility  operating  charges  and  temporary 
facility  charges.  All  costs  are  assumed  to  be  time- stationary. 

Permanent  expansion  costs  will  be  denoted  by  g(-),  where  g(X) 
is  the  cost  of  a permanent  expansion  of  size  X+l.  Notice  that  g(*) 
is  presumed  to  be  a function  only  of  the  expansion  size,  independent  of 
the  level  of  existing  permanent  capacity  and  the  value  of  the  temporary 
facilities  usage  limit,  K.  This  restriction  is  relaxed  in  the  revised 
model  of  [18]. 
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It  Is  important  to  note  that  g(0)  represents  the  cost  of  a unit 
expansion.  Thus,  g(0)  is  comprised  of  the  fixed  (i.e.,  "setup")  expan- 
sion cost  plus  the  marginal  cost  for  the  first  unit  of  increased  capacity. 
Hence,  any  presumption  regarding  a specific  functional  form  (e.g.,  con- 
vexity) for  g over  [0,oo)  will  not  preclude  the  existence  of  a fixed 
expansion  cost. 

With  regard  to  other  costs,  the  temporary  facility  charges  are  of 
principal  interest  in  this  study.  In  order  to  provide  sufficient  generality 
with  regard  to  these  costs,  the  permanent  facility  operating  charges 
will  be  presumed  to  be  proportional  to  the  amount  of  permanent  capacity  used. 
The  following  theorem  indicates  the  modelling  simplifications  that  result 
from  this  assumption. 

Theorem  1.1.  Suppose  that  permanent  facility  operating  charges  are  pro- 
oortional  to  the  amount  of  permanent  capacity  used,  at  rate  y per  unit- 
time. Then  an  equivalent  cost  model  is  given  as  follows: 

( i)  When  permanent  facilities  satisfy  all  demand,  no  operating  charges 
are  incurred. 

(ii)  When  temporary  facilities  are  necessary  (i.e.,  k < 0)  marginal 

costs  equal  to  the  temporary  facility  charges,  less  an  adjustment 
y|k|  = -yk  > 0 per  unit-time,  are  incurred. 

Proof.  Consider  an  arbitrary  demand  pattern  and  an  arbitrary  expansion 
policy.  Let  £( t)  denote  the  demand  at  time  t and  let  k(  t)  denote 
the  spares  level  at  time  t,  t > 0.  Let  ^(  • ) denote  the  0-1  nonnegativity 
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indicator  function:  ^(k)  =0  if  k < 0 and  -^(k)  = 1 if  k > 0. 

Let  x(')  denote  the  complement  of  x('):  x(  • ) = 1 - x(')*  Finally, 
let  g denote  all  future  discounted  expansion  costs  and  temporary  facility 
charges.  Then  the  total  discounted  cost  C is  given  by 


C = 3+/  {x(k(  t))  fi(  t)  r + X(k(  t))  [*>(  t)- |k(  t)  | ]y)  e"rt  dt  (1.1) 

0 


= 3 + Y f £>(  t)  e-rt  dt  - y / xik(  t))  |k(  t)  |e"rt  dt  . ( 1.2) 

0 0 


The  first  term  of  the  integral  in  (1.1)  represents  the  permanent  facility 
operating  charges  incurred  when  the  spares  level  (k(t))  is  nonnegative; 
the  second  term  represents  the  permanent  facility  operating  charges 
incurred  when  -k(t)  >0  customers  must  be  served  by  temporary  facilities. 

The  first  integral  of  (1.2)  is  constant  upon  taking  expectations 
over  the  probability  distribution  of  J9( • ) (independent  of  the  expansion 
policy  chosen).  The  second  integral  of  (1.2)  can  be  accounted  for  in  the 
expected  discounted  cost  minimization  by  subtracting  an  adjustment 
y|k|  = -yk  > 0 per  unit-time  whenever  k(t)  < 0 (i.e.,  whenever 

x(k(c))  = 1).  □ 


The  temporary  facility  charges  will  be  assumed  to  be  dependent  on 
the  value  of  the  spares  level  k and  the  value  of  the  temporary  facilities 
usage  limit  K,  but  independent  of  the  level  of  existing  permanent  capacity. 
Given  these  cost  assumptions.  Theorem  1.1  will  allow  the  construction  of 
Models  I and  II  to  proceed  independent  of  the  actual  level  of  existing 
permanent  capacity. 
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The  presumed  independence  of  the  permanent  operating  costs  from 
the  permanent  facility  capacity  may  not  be  that  restrictive,  since  pro- 
portional fixed  operating  charges  which  depend  solely  on  the  permanent 
capacity  level  can  be  included  in  the  expansion  cost  function  g.  To 
illustrate,  suppose  that  the  actual  permanent  operating  costs  include 
a proportional  fixed  charge  of  v|/  per  unit- time  for  each  unit  of  per- 
manent capacity  available.  Let  Y^  denote  the  initial  permanent  capacity. 

The  expected  discounted  cost  resulting  from  the  t charge  on  the  initial 

00 

* IT  t - ]_ 

capacity  is  then  / YQ\|r  e dt  = \|rr  YQ;  this  cost  is  independent  of 
the  expansion  policy  chosen.  The  expected  discounted  costs  resulting 
from  the  ij/  charge  on  future  capacity  expansion  increments  can  then  be 
accounted  for  by  defining  g as 

00 

g(X)  = g (X)  + / (x+l)t  e'rt  dt 
0 

* ge(X)  + \|ir”  1(X+1)  , 

where  gg  represents  the  actual  expansion  function  alone. 

All  restrictions  (except  time-stationarity)  on  the  facility 
charges,  both  permanent  and  temporary,  are  also  relaxed  in  the  revised 
model  of  [18]. 
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1.6.  Related  Results  and  the  Inventory  Analogy 

In  the  area  of  capacity  expansion,  finite-horizon  models  for 
deterministic  demand  growth  are  most  prevalent.  The  models  of  Manne 
and  Veinott  [13]  and  Erlenkotter  [ 3]  are  representative  of  deterministic 
approaches  employing  dynamic  programming  techniques.  The  models  of 
Bergendahl  [ 1 ] and  Rogers  [ 14 ] are  representative  of  deterministic 
approaches  using  mathematical  programming  solution  techniques.  [ 12 ] 
is  a comprehensive  reference  on  several  deterministic  models. 

The  study  of  Manne  [ 11 ] constitutes  the  first  capacity  expansion 
model  recognizing  stochastic  demand.  In  the  Manne  model,  demand  growth 
is  assumed  to  behave  according  to  a Weiner  diffusion  process.  Expansion 
costs  are  given  by  a concave  power  function: 

g(X)  = aXb,  a > 0,  0 < b <1.  (1.3) 

The  model  alsc  permits  backlogs  at  a proportional  penalty  rate.  Utilizing 
the  Laplace  transform  of  the  presumed  demand  distribution,  Manne  derives 
an  integral  equation  in  X and  K.  This  equation  is  then  solved 
numerically  for  all  feasible  pairs  (X,K)  in  order  to  obtain  minimum 
expected  discounted  cost  parameters  (X*,K*). 

More  recently,  stochastic  capacity  expansion  models  for  telephone 
cable  and  pair-gain  analysis  (see  Example  1.4  of  Section  1.3)  have  been 
reported  in  related  papers  by  Freidenfelds  [ 53,  [ 6 ],  and  Koontz  and 
Shipley  [l0].  Freidenfelds  utilizes  the  Poisson  demand  model  used  here. 
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The  Freidenfelds  model  corresponds  to  Example  2.5  of  Chapter  2 for  the 
special  cases  of  K = 0 and  K = L when  the  expansion  cost  function  is 
given  by 


g(X)  = aX  + b . (1.4) 

By  considering  the  Laplace  transforms  of  first-event  times  (see  Section  2. 
Freidenfelds  derives  a functional  in  X for  the  expected  total  dis- 
counted cost.  This  functional  is  then  minimized  by  function  iteration 
(see  Section  4.4).  Freidenfelds'  results  were  the  stimuli  for  the 
author's  own  investigations,  the  preliminary  findings  of  which  were 
first  reported  in  [10]. 

In  [ 10 ] ^ the  author  presents  a model  corresponding  to  Example  2.5 
of  Chapter  2 when  the  expansion  cost  function  is  given  by  (1.4).  Also 
in  [10],  Koontz  reports  steady-state  results  for  the  same  problem  when 
demand  is  generalized  to  a birth-and-death  process  and  the  minimization 
objective  is  average  cost  per  unit-time. 

Other  related  models  can  be  found  in  the  area  of  single- item 
inventory  theory.  Manne  [ll ] was  the  first  to  note  the  similarity  between 
the  recurrent  (X,K)  expansion  policies,  as  treated  here,  and  (S,s) 
inventory  policies.  Figure  1.4  illustrates  the  near-analogy.  The 
spares  level  (k)  corresponds  to  inventory  on-hand;  a nonnegative  spares 
level  corresponds  to  an  absence  of  inventory  backlogs  and  a negative 
spares  level  corresponds  to  the  existence  of  backlogs.  The  temporary 
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facilities  usage  limit  (K)  corresponds  to  an  inventory  backlog  limit 
(i.e.,  K+l=-s).  The  arrival  process  in  the  expansion  model  corresponds 
to  item  demands  in  the  inventory  model,  while  the  departure  process 
corresponds  to  item  returns.  The  capacity  expansion  function  is  analogous 
to  the  inventory  ordering  (or  production)  cost  function  (i.e.,  X+l  = S-s) 
and  the  temporary  facility  charges  correspond  to  inventory  backlog 
penalties. 

Unfortunately,  the  analogy  is  deficient  with  regard  to  permanent 
facility  operating  charges  and  inventory  holding  charges.  In  the  absence 
of  intervening  expansion,  the  permanent  facility  operating  charges  should 
realistically  decrease  as  the  spares  level  (k)  increases,  since  fewer 
permanent  facilities  are  used.  But,  in  the  absence  of  intervening  orders, 
the  inventory  holding  charges  should  realistically  increase  as  the  on-hand 
inventory  ( k)  increases.  Thus,  the  realistic  behavior  of  permanent 
facility  operating  charges  and  inventory  holding  charges  are  precisely 
opposite,  and  the  analogy  fails.  However,  if  the  permanent  facility 
operating  charges  are  constant  with  regard  to  k,  then  the  analogy  follows 
for  single-item  inventory  models  (with  returns)  having  constant  inventory 
holding  costs,  independent  of  the  level  of  on-hand  inventory.  Referring  to 
Theorem  1.1,  it  follows  that  the  inventory  analogy  also  holds  (with 
inventory  holding  costs  zero)  whenever  the  permanent  facility  operating 
charges  are  proportional  to  the  permanent  capacity  utilized,  provided 
that  the  proper  adjustments  are  made  to  the  temporary  facility  charges 
(i.e.,  the  backlog  penalties).  To  illustrate,  suppose  that  the  permanent 
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facility  charges  are  proportional  at  rate  y per  unit-time  and  that  the 
temporary  facility  charges  are  also  proportional  at  rate  y ' per  unit- 
time, for  each  unit  of  temporary  capacity  utilized.  Then,  employing 
Theorem  1.1,  the  correct  inventory  analogy  is  inventory  holding  costs 
zero  and  backlog  penalty  equal  to  (y'-y)  per  unit-time,  for  each  item 
backlogged.  Thus,  under  the  hypothesis  of  Theorem  1.1,  the  analogy  holds 
for  single-item  inventory  models  with  returns  and  backlogs,  provided 
that  the  inventory  holding  costs  are  set  to  zero  and  the  backlog  penalty 
is  properly  assigned. 

Although  the  literature  for  single-item  inventory  models  is  quite 
extensive  (see  [20]),  the  literature  for  single-item  inventory  models  with 
returns  is  not.  In  fact,  the  literature  for  single- item  inventory  models 
with  both  returns  and  backlogs  is  apparently  null.  Stochastic  inventory 
models  permitting  returns  include  those  of  Tainiter  [19],  Whisler  [21 ], 
and  Heyman  and  Hoad  ley  [7],  [8],  Although  the  demand  processes  of  these 
models  are  similar  to  the  Poisson  model  employed  here,  backlogs  are  not 
permitted.  Furthermore,  all  costs  are  assumed  proportional,  with  the 
exception  of  [21]  where  convex  ordering  functions  are  permitted  over 
discrete  time  periods.  Because  of  the  simplifications  promoted  by 
proportional  costs  and  the  absence  of  backlogs,  the  techniques  employed 
to  analyse  these  models  are  quite  different  than  those  utilized  here. 
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1.7 • The  Untruncated  Demand  Assumption 

As  noted  earlier,  the  demand  process  is  characterized  as  the 
difference  between  two  simple  independent  Poisson  processes.  As  with 
any  mathematical  model,  the  practitioner  must  proceed  with  caution  prior 
to  implementing  the  models  introduced  here  in  an  actual  application. 

However,  there  is  one  assumption  made  here  regarding  the  demand  process 
that  will  always  differ  from  the  actual  case.  Namely,  the  demand  process 
is  not  entirely  legitimate  because  it  allows  negative  demands.  That  is, 
depending  on  the  initial  level  of  demand  $(0)  and  the  arrival  and 
departure  rates,  there  is  some  probability  that  customers  will  depart 
when  there  are  no  customers  in  the  system^  This  anomaly  of  the  process 

j 

will  be  referred  to  as  the  untruncated  assumption,  reflecting  the  fact 
that  the  models  do  not  truncate  departures  when  demand  becomes  zero. 

The  purpose  of  this  section  is  to  demonstrate  that  the  untruncated  assump- 
tion will  not  significantly  affect  the  derived  optimal  operating  parameters 
in  most  cases  and  to  identify  the  few  cases  where  it  might  adversely 
distort  results.  With  regard  to  these  later  cases,  the  practitioner  might 
well  consider  a revised  model  on  the  order  of  that  introduced  in  [5], 
which  provides  truncation. 

The  reason  for  making  the  untruncated  assumption  is  one  of 
mathematical  simplicity.  Coupled  with  other  basic  assumptions,  the 

untruncated  assumption  will  allow  model  constructions  that  need  not 
explicitly  account  for  the  total  level  of  permanent  capacity  available 
and  the  total  system  demand  level.  Mathematically,  the  models  can  then 
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be  characterized  by  a single  set  of  equations,  as  opposed  to  a separate 
set  of  equations  for  each  possible  level  of  permanent  capacity. 


I 

! 


The  models  of  Manne  [ll]  and  Freidenfelds  [6]  also  possess  this 
anomaly  (the  anomaly  is  not  present  in  the  inventory  analogy) . An  argument 
was  presented  in  [ll]  that  purported  to  demonstrate  that  the  untruncated 
assumption  did  not  affect  the  optimal  operating  parameters  at  all. 

Briefly,  the  argument  can  be  paraphrased  as  follows:  "since  no  costs 

are  incurred  when  demand  is  negative,  the  total  costs  of  the  untruncated 
model  will  be  identical  to  those  of  the  actual  truncated  case.  Therefore, 
the  operating  parameters  derived  from  the  untruncated  model  will  be  optimal 
for  the  actual  truncated  case."  Concerning  the  case  at  hand,  observe 
that  the  spares  level  (k)  will  be  positive  whenever  demand  becomes 
negative.  Hence,  in  view  of  Theorem  1.1  it  can  be  assumed  that  no  costs 
are  incurred  whenever  demand  becomes  negative  in  the  models  considered 
here.  Thus,  the  Manne  argument,  if  valid,  would  also  hold  for  the  model 
at  hand. 

Unfortunately,  the  argument  put  forth  by  Manne  has  two  flaws. 
Firstly,  once  demand  becomes  negative,  there  is  some  probability  that 
it  will  remain  so  thereafter;  for  these  sample  paths,  the  total  cost 
in  the  untruncated  model  can  differ  from  the  actual  truncated  situation. 
Secondly,  even  though  total  costs  may  be  the  same  for  many  of  these 
sample  paths,  total  discounted  costs  ( the  criteria  used  both  here  and 
in  [11  ])  may  not  be.  In  fact,  for  those  sample  paths  where  demand 
becomes  negative,  future  expenditures  occur  later  (if  at  all)  than  in 
the  actual  truncated  case.  Hence,  the  total  discounted  cost  for  these 
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estimate  of  the  actual  situation.  It  will  now  be  demonstrated  that  this 
estimate  is  quite  adequate  for  most  cases. 

As  shown  in  [ 6 ],  the  expectation  and  variance  of  demand  £(•) 
are  given  by 


E[£(t)]  = fl(0)  + At(  p-  q)  , t > 0 


Var[s( t) ] = At  , t > 0 


(1.5) 


(1.6) 


Thus,  when  p > q,  the  mean  demand  level  increases  with  time;  this  will 
be  referred  to  as  the  growth  situation.  When  p < q,  the  mean  demand 
level  decreases  with  time;  this  will  be  referred  to  as  the  negative 
growth  situation.  The  case  p = q will  be  referred  to  as  the  stable 
situation.  The  events  of  interest,  with  regard  to  the  untruncated 
assumption,  are 


e(  t)  = {jfl(  s)  <0  for  some  s S [0, t]},  t>0  . 


First  consider  P{ £(<»)),  the  probability  of  demand  eventually  becoming 
negative.  Given  the  independent  characterization  (A, p)  introduced  in 
Section  1.4  for  the  demand  process,  it  should  be  apparent  that  this  probability 
is  equal  to  the  probability  that  a random  walk  (p,q)  starting  at  the 
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origin  ever  reaches  the  position  -(£(0)  + 1).  From  [ 4 j,  this  probability 
is 


P(  &H) 


(q/p)^0^1  , 

1 , 


p > q 

p < q 


That  is,  the  probability  of  demand  eventually  becoming  negative  is 
( q/p)^°)  + 1 for  the  growth  situation;  for  other  cases,  the  event  is 
certain.  Thus,  if  $(0)  is  sufficiently  large  in  the  growth  situation, 
then  the  proportion  of  sample  paths  exhibiting  negative  demands  will  be 
small  and  the  untruncated  assumption  should  be  viable.  In  fact,  since 
no  costs  are  incurred  while  the  spares  level  is  positive  (Theorem  1.1)  it 
suffices  to  assume  that  $(0)  equals  the  total  permanent  capacity 
initially  available.  Thus,  if  the  initial  permanent  capacity  is  suffi- 
ciently large  in  the  growth  situation,  then  the  untruncated  assumption 
should  be  viable.  To  illustrate,  for  p = .52  (q  = .14-8)  and  £)(0)  = 100, 
P(£(°°))  = 3.5  X 10  \ Therefore,  the  untruncated  assumption  does  not 
appear  unreasonable  for  most  growth  situations  where  there  is  some 
initial  permanent  capacity  to  begin  with. 

Although  P{&(oo)j  = 1 for  non-growth  cases,  the  untruncated 
assumption  is  still  reasonable  for  many  of  these  situations.  This  fact 


results  from  two  basic  observations: 

( i)  if  the  discount  rate  is  not  abnormally  small,  then  a large  portion 
of  the  total  expected  discounted  cost  occurs  over  the  finite 
interval  [0, t];  and 
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(ii)  of  those  sample  paths  contributing  to  the  event  g( t) , only  a portion 
of  these  paths  will  exhibit  discounted  costs  which  differ  from 
their  truncated  counterparts  over  [0, t]. 

Observation  ( i)  follows  from  the  fact  that  the  models  are  time- 
stationary. Since  the  costs  and  the  demand  process  do  not  change  over 
time,  the  expected  expenditures  over  each  time  interval  of  length  t 
will  be  approximately  equal,  for  t sufficiently  large.  Let  A denote 

the  expected  discounted  costs  over  (0  t],  Then  the  total  expected  dis- 

* IT  t » 

counted  costs  over  [O,oo)  can  be  approximated  by  A( 1 + e + e +•••) 
■■  r t — 1 

= A( 1 - e ) . Thus,  the  expected  discounted  costs  incurred  over 

IT  t 

[0,t]  will  account  for  approximately  ( 1 - e"  ) x 100$  of  the  total 

expected  discounted  costs.  To  illustrate,  for  t = 10  and  r = . 15, 

(1  - e"rt)  x 100 $ = 78$;  for  t = 20  and  r = .10,  (1  - e”rt)  x 100$  = 86$. 

Concerning  observation  ( ii) , let 

t = min(  s : f)(  s)  < 0} 

Then  £(t)  = (t  < t}.  As  noted  earlier,  it  suffices  to  assume  that 
fl(0)  equals  the  initial  permanent  capacity  level.  At  time  T,  the 
spares  level  in  the  actual  untruncated  case  is  equal  to  the  permanent 
capacity  level;  a lower  bound  on  the  permanent  capacity  level  at  time  r 
is  jB(0).  By  Theorem  1.1,  no  additional  expenditures  will  occur  until 
the  spares  level  reaches  zero;  a necessary  (but  by  no  means  sufficient) 
condition  for  this  to  happen  is  that  at  least  &(0)  additional  arrivals 
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occur  past  T.  Hence,  of  those  sample  paths  contributing  to  the  event 
£.(  t) , only  those  exhibiting  at  least  &(0)  additional  arrivals  during 
the  interval  (r,  t]  can  possibly  have  discounted  costs  which  differ 
from  the  actual  truncated  case  during  the  interval  [0, t].  Let  Na(s) 
denote  the  total  number  of  arrivals  during  the  interval  [0,s].  Then 
an  upper  bound  on  the  proportion  of  sample  paths  exhibiting  discounted 
costs  which  differ  from  the  actual  truncated  situation  over  [0,  t]  is 
given  by  U(  t)  = P{t  < t,  Nfl(  t)  - Nfl(  r)  >£(0)}. 

Since  the  arrival  process  is  Poisson  at  rate  Xp, 


P{N  (t)  - N (t)  >fi(0)|r  = s < t)  = Z .{AP(t-s)}n  e->,p(  t-s)  _ 

" n=s9(0) 

Therefore,  letting  f and  F respectively  denote  the  probability 
density  and  cumulative  distribution  functions  of  t 


U(t  )=/  f(s)  Z 

0 n=jB(0) 


< F(t)  = p{  e(  t) ) . 


I r-s)  1 ' e"^P(  t-s) 

n!  6 


(1-7) 


No  attempt  will  be  made  here  to  actually  compute  U( • ) . Instead,  F( • ) 
will  be  estimated  using  a normal  approximation.  The  significance  of 
(1.7)  is  that  the  approximation  for  F(.)  should  be  entirely  adequate 

since,  for  J5(0)  sufficiently  large,  U(«)  will  be  many  orders  of 
magnitude  smaller  than  F(>)  (recall  also  that  U(>)  itself  is  a 
rather  loose  upper  bound  on  the  actual  relevant  probabilities). 


ET  .r-- ; 


a 


The  following  lemma  discloses  the  form  of  F and  provides  further 
justification  for  the  validity  of  the  normal  approximation  to  be  used. 


Lemma  1.2.  Let  m = jfl(O)  + 1.  Then 


F<t)  = l ?((J"w2)  P(m_j)/2  q(m+j)/2  l ^e-At,  t >0(1.8) 


j=m  j Mm+j)/2; 


(Note:  ^(m+j)/2^  = ® if  the  parities  of  m and  j differ.) 

Proof:  Let  N( t)  denote  the  number  of  events  (arrivals  and  departues) 

during  the  time  interval  [0,t].  Since  (N( t) , t > 0)  is  Poisson, 


P{N(t)  = n)  = e"At  , 


t > 0 


(1.9) 


Given  N( t)  = n,  we  seek  the  probability  that  the  demand  becomes 
negative  during  the  first  n events.  Let  P[£..)  = P(demand  first 
becomes  negative  at  jth  event).  Obviously,  P(t^)  = 0 f°r  j < m, 
since  at  least  m departures  are  required.  As  shown  in  [ 4 ] (pages 

351-352), 


pip  l _ J2  / m . \ D(m-j)/2  a(m-j)/2 
tejJ  ' j (nH-j)/2  P q 


j > m 


(1.10) 


The  derivation  of  (1.10)  follows  from  the  fact  that  in  order  for  demand 
to  first  become  negative  at  the  jth  event,  there  must  be  precisely 
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(m-j)/2  arrivals,  (m+j)/2  departures,  and  the  demand  must  be  nonnegative 
at  all  previous  event  epochs. 


1 


1 


Thus,  from  ( 1.9)  and  (1.10), 


F(t)  = P{£(t)) 


= Z ( Z P{£j)  P(N(t)  = n] 
n=m  j=m  J 


= Z P(&.)  Z P(N(t)  = n) 
j=m  J n= j 

00 

= Z P{£.)  P{N(t)  > j} 
j=m  J 


_ V ® ( m . ) p(®-j)/2  Q(nn-j)/2  y CKt)n  -At 

- /:mj  ((-j)/2)  p q n: 


Notice  that  ( 1.8)  is  in  essence  of  partial  convolution  of  binomial 
(almost)  and  Poisson  probabilities.  Therefore,  for  m (i.e.,  fi(0)) 
sufficiently  large.  At  > 10,  and  neither  p nor  q close  to  zero,  a 
normal  approximation  to  (1.8)  should  be  applicable  [ 2 ]. 

In  order  to  implement  a normal  approximation,  the  mean  and  variance 
of  t are  required.  For  p = q the  mean  of  t is  infinite.  However, 
it  is  finite  for  the  negative  growth  situation. 


Then,  T=  T-(jg(o)  + l) ' Let  ^ 0 denote  the  probability  density  function 
of  x^.  Then,  by  conditioning  on  the  time  ( s)  and  the  type  (arrival 
or  departure)  of  the  first  event,  the  following  integral  difference 
equation  can  be  written  for  f ( t) : 


Let  h( i)  denote  the  Laplace  transform  of  using  (I.I3), 


h(i)  = E[e"“i]  = / f ( t)  e-rt  dt 

n * 


00  t 

f 0-rt  r ^-As 


= / e / Ae‘  S[pfi_l(t-s)  + qfi+1(t-s)}ds 


= A / e-AS  / 
0 s 


As 

{ 6 t(pfi-l(t"8)  + qfi+1(t-s)}dt  ds 


" * 0 ' **  8 " { e"r(t“)tP«1.1(‘-*)  * <!£1+1( t-1 

00 

= A / e A+r^S{ph(i-l)  + qh(  i+1)  }ds 


= A(  A+r)  {ph(i-l)  + qh(  i+1) ) 


1 <-l. 


where  h(0)  = 1. 

The  solution  to  difference  equations  (1.14)  is  given  in  [ 4 ] 
and  [ 9 ] as 


Using  (1.15) 


F.  r t i - W-itlP-s)  , i~1 l — i 

Mp-q)5  * (P-q)  ^ (p-q)5  A2(p-q)^ 


Finally, 


Var[T.]  = E[x2]  . (E[Ti])2  = — i , 

A (p-q)5 

where  verifies  (1.12)  upon  substitution  of  i = _(jq(o)  + 1). 


Using  the  above  theorem,  the  normal  approximation  for  F( t)  is 


given  by 


/t  . iXQM  x 

f^A(q-p)  ^ 

E/ 


= F(t) 


A2(  q-p)5 


At(  q-c 


f£(o)±i) 


( jfi(  0)  + l) 


_l-2p)3/2  - (l-2p)l/2  (afoUtt 


(i(o)+l)' 


, q > p 


, P < .5  , (1.16) 


where  $ denotes  the  cumulation  standard  normal  distribution  function 
(also  recall  that  q = 1-p) . Using  standard  normal  tables,  the  right 
hand  side  of  (1.16)  will  not  exceed  .01  whenever  the  argument  does  not 


exceed  -2.326.  That  is,  whenever 
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(1.17) 


Xt(l-2p)5/2  - (l-2p)1/2  (fi(O)i-l)  + 2.326  (jB(0)+1)  1/2  < 0 . 


Using  the  quadratic  formula,  ( 1. I7)  will  be  satisfied  whenever 

*(o)  > i + '2'?26  + jfUp1)  * 4 ■t('I'2pl2j  ■ 


(1.18) 


It  is  a simple  matter  to  see  that  (1.18)  increases  with  X and  t , and 
decreases  with  p.  Table  1.1  is  a small  tabulation  of  (1.18)  for 
X = 100,  .35  < p < .45,  and  t = 10  and  20. 


t 


X = 100 

10 

20 

.55 

385 

715 

.40 

290 

520 

.45 

1 

207 

336 

Table  1.1:  Lower  Bounds  on  &(0) 


for  F(t)  < .01 


To  illustrate,  for  p > .4,  Table  1.1  indicates  that  if  £(0)  >520 
and  X < 100,  then  less  than  1$  of  the  untruncated  sample  points  will 
exhibit  negative  demands  over  the  first  20  years  (and  approximately  Q6t 
of  the  total  discounted  costs  are  determined  in  the  first  20  years  if 
r > . 1) . To  see  that  Table  1.1  is  not  unreasonable,  consider  the 
statistic  | = X/f)( 0).  | represents  the  amount  of  expected  "churning" 

or  "turn-over"  in  the  system  relative  to  the  initial  demand  jfi(O).  That 
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is,  | approximates  the  proportion  of  initial  demand  involved  with 
transactions  (arrivals  and  departures)  during  a single  year.  Thus, 
jfl(O)  > 520  for  A = 100  requires  that  expected  annual  turnover  must  be 
somewhat  less  than  20%  in  order  to  insure  that  F(20)  < .01. 

Therefore,  in  nongrowth  situations,  the  untruncated  assumption 
appears  to  be  viable  for  situations  where  £(0)  is  sufficiently  large 
relative  to  the  annual  event  rate  A;  that  is,  when  annual  turnover  is 
sufficiently  small.  This  must  be  the  case  in  most  nongrowth  situations 
where  future  permanent  expansion  is  a realistic  possibility.  For,  if 
A is  large  relative  to  j9(0)  and  p < q,  then  E[fi(t)]  tends  toward 
0 rather  rapidly  according  to  ( 1.5)  i in  this  case,  future  capacity 
contraction  --  not  expansion  --  is  the  alternative  requiring  investiga- 
tion. A similar  remark  holds  for  cases  where  p is  small  (e.g.,  P < .35) 

In  conclusion,  the  untruncated  demand  assumption  appears  to  be 
viable  in  most  cases  where  $(0)  is  sufficiently  large  and  future 
capacity  expansion  is  a realistic  possibility.  For  growth  situations 
(p  > q) , the  assumption  will  be  viable  for  nearly  all  applications  with 
$(0)  >0  (i.e.,  existing  initial  permanent  capacity).  For  nongrowth 

situations,  the  assumption  appears  to  be  viable  for  applications  where 
additional  permanent  expansion  is  a realistic  future  possibility. 


Li> 


..8.  Notation 

As  introduced  previously  in  this  chapter,  $(  t)  will  denote  the 
demand  at  time  t > 0,  A will  denote  the  event  rate  with  arrival  prob- 
ability p and  departure  probability  q,  k will  denote  the  spares 
level,  K will  denote  the  temporary  facilities  usage  limit,  X+l  will 
denote  the  expansion  size  and  g will  denote  the  permanent  expansion 
cost  function,  where  g(X)  represents  the  cost  for  an  expansion  of 
size  X+l. 

IR  will  denote  the  space  of  all  real  vectors  having  m components, 
m >_  1.  Only  row  vectors  will  be  utilized.  Given  U,V  £ ]Rm,  the  dot 
product  U*V  will  be  used  to  denote  vector  multiplication.  RmXn  will 
denote  the  space  of  all  real  matrices  having  m rows  and  n columns. 
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CHAPTER  2 


MODEL  I:  THE  NON-MODULAR  CASE 

2. 1.  Introduction 

This  chapter  characterizes  the  untruncated  constant  Poisson  model 
under  two  basic  assumptions: 

Assumption  2. 1.  Permanent  facility  operating  charges  are  proportional  to 
the  amount  of  permanent  capacity  used,  at  rate  per  unit- time. 

Assumption  2.2.  Permanent  facilities  are  always  used  in  preference  to 
temporary  facilities  and  temporary  facility  charges  are  incurred  only 
when  there  is  excess  demand  to  be  served  (i.e.,  k < 0) . Furthermore, 
these  charges  are  functions  solely  of  the  excess  demand  (-k  > 0),  for 
a given  limit  (K)  on  temporary  facility  usage. 

Assumption  2. 1 is  a restatement  of  the  hypothesis  for  Theorem  1. 1. 
Therefore,  using  Theorem  1.1,  the  construction  of  Model  I will  proceed 
under  the  equivalent  assumption: 

Assumption  2,1'. 

( i)  When  permanent  facilities  satisfy  all  demand,  no  operating  costs 
are  incurred. 

( ii)  When  temporary  facilities  are  necessary  (i.e.,  k < 0)  marginal  costs 
equal  to  the  temporary  facility  charges,  less  an  adjustment 
T ^ I k | = -y^k  > 0 per  unit-time,  are  incurred. 
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Assumption  2.2  presumes  a known  unique  correspondence  between 
temporary  facility  costs  and  the  value  of  the  spares  level  k < 0,  given 
the  usage  limit  K > 0.  It  is  demonstrated  in  [17]  that  this 
assumption  does  not^  in  general^  hold  for  the  case  of  modular  temporary 
facilities.  However^  the  assumption  appears  to  be  viable  for  any  situa- 
tion where  temporary  facility  costs  are  not  a function  of  the  past  temporary 
facility  usage  pattern.  For  instance^  referring  to  the  inventory  analogy 
of  Chapter  1}  Assumption  2.2  is  always  presumed  with  regard  to  backlogging 
charges . 

These  assumptions  will  be  treated  more  precisely  in  the  next 
section.  In  the  meantime^  the  following  examples  are  offered  to  illustrate 
situations  where  they  hold. 


EXAMPLES 

Example  2. 1 ; General  Single  Function.  Suppose  that^  given  the  value  K, 
temporary  facility  charges  are  given  by  a function 

w^(k^t)  = discounted  temporary  facility  costs  for  serving 
| k | customers  over  a length  of  time  t > 0, 
k < 0. 


In  this  case;  the  marginal  costs  can  similarly  be  written  as 


costs  given  by  any  of  the  previous  examples>  an  instantaneous  charge 
£ 

is  incurred  whenever  an  additional  unit  of  temporary  capacity  is 
used  (that  is^  each  time  that  k decreases  whenever  k < 0).  This  case 
differs  from  the  general  modular  case  because  the  instantaneous  charge 
is  solely  a function  of  k and  is  independent  of  the  past  temporary 
facility  usage  pattern. 
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Example  2.j:  Restricted  Modular  Case.  As  noted  previously,  Assumption  2.2 

does  not  hold  in  the  general  modular  case.  However,  the  assumption  will 
hold  _if  the  utilization  of  temporary  modules  is  required  to  track  the 
demand  pattern.  Recalling  that  each  module  provides  L units  of  temporary 
capacity,  there  is  a one-to-one  correspondence  between  the  spares  level 
(k)  and  the  number  of  modules  used  (n)  under  this  restriction: 

( 

0 , k > 0 

n = / 

min{i  : iL  > -k)  = [ , k < 0 

\ 

where  [•]  denotes  integer  truncation. 

This  situation  can  be  represented  in  a transition  diagram,  as 
shown  in  Figure  2.1.  In  the  diagram,  spares  levels  are  depicted  vertically 
and  spares  levels  using  the  same  number  of  temporary  modules  are  aligned 
in  columns.  Each  circled  number  denotes  the  number  of  temporary  modules 
necessary  for  the  corresponding  spares  level  given  horizontal  to  the  circle 
in  the  right-hand  margin.  Transitions  occurring  due  to  departures  are 
denoted  by  solid  lines,  while  transitions  occurring  due  to  arrivals  are 
denoted  by  dashed  lines. 

Connection  and  disconnection  (i.e.,  installation  and  removal) 
charges  are  denoted  by  triangles.  Whenever  all  modules  are  fully  used 
(k  = -nL  for  some  integer  n)  and  an  arrival  next  occurs,  two  possibilities 
occur.  If  k = -K  (=  -NL  for  some  integer  N) , then  a permanent  capacity 
expansion  is  undertaken  and  the  N modules  are  disposed  of;  otherwise 


(2.1) 
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(n  < N),  a new  module  is  connected  at  cost  c,  increasing  the  total 
number  of  modules  in  use  to  n+1.  Similarly,  whenever  a module  is  being 
used  to  serve  a single  customer  (i.e.,  k = -(n-l)L-l  for  some  n > 1)  and 
a departure  next  occurs,  a module  is  disconnected  at  cost  d,  decreasing 
the  number  of  modules  in  use  to  n-1. 

In  addition  to  the  instantaneous  connect  and  disconnect  charges,  a 
cost  rate  of  jt  per  unit- time  is  incurred  for  each  module  in  use.  Given 
the  correspondence  (2.1),  adjusted  cost  rate  [L  ^(-k+L-1)]  (n  + T-jk) 
per  unit- time  is  incurred  whenever  k < 0. 

Given  "free  rein"  on  temporary  module  utilization,  it  is 
shown  in  [17]  that  optimal  disconnections  will  usually  not  track 
the  demand  pattern.  Therefore,  Model  I is  only  applicable  for  the  modular 
case  when  the  restriction  of  module  use  to  the  demand  pattern  is  an 
imposed  operating  constraint.  This  situation  can  arise  when  temporary 
modules  are  dispatched  from  a central  facility  which  services  a number  of 
systems,  including  the  one  presently  being  modeled.  If  the  overall  supply 
of  modules  is  limited,  then  the  restriction  of  module  use  to  the  demand 
pattern  may  be  an  imposition  dictated  by  the  central  facility. 

The  restriction  can  also  arise  when  the  retention  of  unnecessary 
modules  is  aesthetically  unpleasing  to  some  individual  or  interest  group. 
For  example,  suppose  that  mobile  ("barracks")  facilities  were  located  on 
a school  campus  to  augment  classroom  capacity.  If  congestion  at  the 
school  temporarily  declines,  then  the  school  board  may  have  difficulty 
convincing  its  constituency  that  retention  of  the  mobile  facilities  is 
optimal. 


1 


b2 


Remark-Disposal  Costs.  Recall  that  whenever  the  temporary  facilities 
usage  limit  is  reached  (i.e.,  k -K)  and  an  arrival  next  occur,  a 
permanent  capacity  expansion  of  size  X+l  is  undertaken  at  cost  g(X), 

X > K.  When  this  happens,  the  spares  level  increases  from  -K  to 
X-K  > 0 and  (by  Assumption  2.2)  permanent  facilities  are  again  used  to 
satisfy  all  demand.  Thus,  when  an  expansion  occurs,  there  may  be  abnormal 
costs  (or  revenues)  for  "turning  off"  or  disposing  of  the  temporary 
facilities.  Thus,  in  all  of  the  above  examples,  there  may  be  an  additional 
cost  (or  revenue)  0g(K)  whenever  k = -K  and  an  arrival  next  occurs. 


2.2.  Transition  Equations  for  the  Expected  Discounted  Costs 

C^(X,K)  will  denote  the  expected  total  discounted  cost  when  the 
initial  spares  level  is  k,  the  temporary  facility  usage  limit  is  K, 
and  the  permanent  capacity  expansion  size  is  X+lj  X > K > 0,  k > -K. 
Similarly,  Fk(K)  will  denote  the  expected  incremental  discounted  costs 
(until  the  next  arrival  or  departure)  while  the  spares  level  has  value  k, 
when  the  temporary  facilities  usage  limit  is  K > 0;  k > -K.  Given  these 
designations.  Assumptions  2.1'  and  2.2  can  now  be  restated: 

(a)  Fr(K)  = 0,  k > 0,  K > 0; 

(b)  F^(K),  0 > k > -K,  K > 0,  are  (uniquely)  known,  and 

(c)  F^K)  includes  an  (expected)  adjustment  term 

/ r.kr'  ( 1 - e~rt)  Ae‘At  dt  = (X+r)‘  y^k  , k < 0 . 
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Recall  that  p is  the  probability  of  arrival  and  q = 1-p  is 


the  probability  of  departure.  Denote: 


OC  = ^q(X+r)  1 and  0 = Ap(A+r)_1 


Definition  2.1.  The  expected  incremental  costs  {F^(K) } will  be  said 
to  be  constant  in  K if  Fk(K)  = FR,  0 > k > -K+l  (independent  of  K) 
for  all  K > 0.  In  this  case,  the  expected  incremental  cost  at  spares 
level  -K  will  be  denoted  F + 0 (K). 

The  ramifications  of  the  above  definition,  revealed  in  Section  2.k, 
are  relatively  minor.  However,  in  actual  practice,  it  appears  likely 
that  many  of  the  examples  in  the  previous  section  will  often  have 
expected  incremental  costs  that  are  constant  in  K.  Given  the  apparent 
prevalence  of  this  property,  the  separate  term  0e(K)  will  be  included 
in  the  model  development  so  that,  without  loss  of  generality,  it  can  be 
assumed  that  F ^(K)  = F R when  the  incremental  costs  are  constant  in 
K.  If  the  incremental  costs  are  not  constant  in  K,  then  0e(')  can 
be  taken  as  identically  zero. 


EXAMPLES  (continued) 


wk(’>')  = w('»")  For  K > 0,  then  the  incremental  costs  are 

constant  in  K. 

Example  2.2.  Conditioning  on  the  time  of  the  next  event  gives 

00 

Fk(K)  = / (r|(k)  + rik)  r_1(l  - e-rt)  Ae‘At  dt 

= (r2(k)  + rLk)  (A+r)"1  , k < 0 . 

£ 

If  Y"2(k)  = y2(k)  for  all  K > 0,  then  the  incremental  costs  are  constant 
in  K. 

K K 

Example  2.3.  When  r2(  k)  = -y-gk,  the  previous  expression  gives 
Fk  = (Tj_"T2)k  (A+r)  f k < 0.  If  = y*g  for  all  K > 0f  then  the 
incremental  costs  are  constant  in  K. 

Conditioning  on  the  time  of  the  next  event  and  the  prob- 
that  the  next  event  is  an  arrival  gives 

00 

I (\(k,t)  + r” 1 r,k(l  - e"rt)  + e"rt)  *e'At  dt 
0 ^ 

00 

/ wR(k,t)  Ae  t dt  + (A+r)-1  y k + Py*  , 


0 > k > -K 


p 


! 


00  -V 

F_r(K)  = / wK(-K,t)  Ae'At  dt  + (A+r)'1  r^k  . 


Note  that  the  instantaneous  charge  y^  is  not  incurred  when  k = -K 

0 * 

since  an  arrival  at  that  point  triggers  an  expansion. 

£ 

If  wlJ  * > = ' > ‘)  and  = ry  for  a11  K > then  denote 


F.  = / w(k,t)  Ae‘At  dt  + (A+r)"1  r,k  + £r,  , 
0 


k < 0 


In  this  case;  F^(K)  = F 0 > k > -K.  Hence,  the  expected  incremental 
costs  are  constant  in  K.  Thus,  the  expected  incremental  cost  when 
k = -K  will  be  denoted  as  F + (30  ( K)  where  the  function  0 (•) 
includes  an  adjustment  -y?. 


Example  2.5.  Using  the  correspondence  (2.1),  the  expected  incremental 
usage  costs  are  given  by 


-k+L-1 


■rt.  , -At 


Rk  = / ( [~  l ' ' ] n + Fjk)  r'  ( 1 - e'  ) Ae  dt 


= « + rLk)  U+r)-1  , 


k < 0 . 


For  connect  and  disconnect  points,  conditioning  on  the  both  the  time 
and  type  of  the  next  event  gives 


. 


Rk> 

Rk+Od, 

if 

if 

mod( -k,L)  ^ 
mod(-k,L)  = 

o,1 

1 

y 

k < 0 

Fk'i 

Rk  + 0C> 

if 

raod(-k,L)  = 

0 

Thus^  Fk(K)  = Fk,  0 > k > -K  (=  - NL) , so  the  expected  incremental  costs 
are  constant  in  K for  the  modular  case.  Hence,  the  expected  incremental 
cost  when  k = -K  will  be  denoted  F ^ f where  the  function 

^e(‘)  includes  an  adjustment  term  -c.  Note  that  0g(K)  likely  also 
includes  a term  Kd/L  to  provide  for  disconnection  of  the  N = K/L 
modules  when  a permanent  expansion  occurs. 

By  conditioning  on  the  time  and  type  of  the  next  event,  the  cost 
transition  equations  are 


Ck(X,K)  = Fk(K)  + f Ae'At(qCk+1(X,K)  + pCk_L(X,K))  e‘rt  dt 
= Fk(K)  + OCk+1(X,K)  + 3Ck_L(X,K) 

“ktlf*^)  + eCk-l(X'K>>  k>0  (2-3) 


= 


|^Fk(K)  + QCk+1(X,K)  + pCk  l(X,K),  0 > k > -K  (2.4) 


Notice  that  no  upper  bound  is  imposed  upon  k in  (2.3).  This  corresponds 
to  the  "untruncated  demand  assumption"  discussed  in  Chapter  1.  Therefore, 
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the  comments  of  Chapter  1 regarding  the  viability  of  this  assumption 
are  applicable  here. 

When  k = -K  and  an  arrival  next  occurs,  a permanent  capacity 
expansion  of  size  X+l  (X  > K)  is  undertaken.  Whenever  this  happens, 
the  spares  level  increases  to  the  value  X-K  > 0.  Therefore,  by  again 
conditioning  on  the  next  event  time  and  type,  the  boundary  transition 
equation  becomes 


C_k(X,k)  = F_k(K)  + / *e‘AC{qC_K+1(X,K) 

+ p(0e(K)  + g(X)  + Cx_k(X,K)))  e"rt  dt 

= F.k(K)  + QC_k+i(X,K)  + P(0e(K)  + g(X)  + Cx_k(X,K)).  (2.5) 

Notice  that,  for  any  scalar  v,  F (K)  + ^v  and  0 (K)  - v may  be 
readily  substituted  for  F_K(K)  and  0 ( K) , respectively,  in  (2.5). 

This  fact  will  be  referred  to  as  the  Equivalent  Formulation  Property. 


2.J.  Solutions  for  the  Transition  Equations,  k > 0 

The  special  nature  of  the  linear  equations  (2.5)  is  well-known 
and  solutions  are  given  in  [ 4 ] and  [ 9 ].  For  the  sake  of  completeness, 
the  appropriate  results  are  repeated  here  in  both  analytical  and  prob- 
abilistic contexts.  The  following  lemma  states  some  prerequisite  rela- 
tionships regarding  the  coefficients  a and  {3. 
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Lemma  2.1.  For  a and  £ given  by  (2.2): 


f. 


! 

;i 


(i)  a > o,  p > o, 

( ii)  a+p  < i, 

( iii)  a*3  < 1/4. 

Proof,  (i)  follows  since  A,  p,  q and  r are  all  positive,  (ii)  follows 
also  from  positivity  since  04-P  = A(A+r)_1.  The  function  p(l-p)  attains 
a global  maximum  of  1/4  when  p = 1/2.  Thus,  Qp  = A (A+r)  pq  < pq 
= p(l-p)  < 1/4,  which  is  (iii).  □ 

The  following  theorem  provides  the  analytical  solutions  to  equations 

(2.3). 

Theorem  2.2. 


Ck(X,K)  = Zk  Cq(X,K), 


k > 0 , 


(2.6) 


where 


z _ 

2a 


(2.6a) 


Furthermore,  Z is  real  and  £ < Z < a+p. 


I 

1 
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Proof.  It  is  easily  verified  that  C^X^K)  = v Cq(X^K),  k > 0,  satisfies 

2 

(2.3)  if  v is  a root  of  the  quadratic  expression  Ov  - v + f3  = 0. 

This  quadratic  has  two  roots: 


1 -1TTT&  , w /TTT5p 

2a  2a 


(2.7) 


By  Lemma  2. 1,  both  roots  are  real  and  positive,  Also, 


OfP  < 1 =>  p < l-a 


=>  Uap  < i+a(  l-a) 

=>  l-Uap  > l-i+cM+a2  = (l-2a)2  = (2a- 1)2 
=»  Vl-4ap  > l-2a  and  vl-i+ap  > 2a- 1 


=*>  1 - Vl-Uap  < 2a  and  1 + /l-J^ap  > 2a 


=>  Z < 1 and  Z ' > 1 . 


To  verify  that  Z is  the  correct  root  for  the  case  at  hand,  we 
rely  on  economic  reasoning.  Specifically,  since  no  charges  are  incurred 
when  k > 0,  no  costs  can  accrue  until  the  first  time  that  the  spares 
level  becomes  nonpositive.  Therefore,  since  the  model  is  time- stationary 
it  must  be  less  costly  to  start  with  a larger  positive  spares  level;  i.e. 
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L a. 


T""»'  • ‘ 


Hence 


Ck(X.K)  < C0(X,K),  k > 0.  Since  Z ' > Cfc(X,K)  ^ Z'k  C0(X,K). 
Ck(X,K)  = Zk  Cq(X,K),  k >0. 

Regarding  the  lower  bound  on  Z}  we  have 


( 1-SQp)  = 1-tap  + 4(Qp)2  =»  1-tap  < ( 1-2QP)2  =*  < i_' 


=>  i - /i-tap  > 2ap  =>  z > p 


The  upper  bound  on  Z follows  as 


Z = QZ  +(3  < 0^.(3  } since  Z < 1 . 


Corollary  2.?.  For  X > K > 0,  X ' > K ’ > 0 and  k > 0, 


VX>K)  < C0(X',K')  <=>  Cfc(X,K)  <C(X',K')  . 


Proof.  The  result  is  a direct  consequence  of  the  theorem  sii 
Z > p > 0. 


The  solution  (2.6)  can  also  be  obtained  through  a simple  probabilistic 


argument. 


Definition  2.2.  Let  the  random  variable  ^ denote  the  first  time  that 


the  net  demand  reaches  i: 


c±  = min(t  >0  : 4(t)  - «(0)  = i},  i = 0,  +1,  +2,  ...  . 


Lemma  2_,  4 . Z = E[e  r^l). 
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f 


-rTti 


Proof.  Let  v = E[e  ].  Since  the  demand  process  is  memory  less, 

Tg  = T'  + T"  where  I'  and  T"  are  independent,  T'  ? and  T"  £ T^. 


Therefore. 


Ete'rT2]  - E[-r(I'>T">]  . E[«-rI']  El.-'1'']  , „2  . 


Conditioning  on  the  first  event  type  and  the  first  event  time  ( t)  yields: 


E[e  rTl! 1st  event  at  t is  on  arrival]  = e”rt 


E[e  ^ I Is t event  at  t is  a departure]  = e~rt  E[e  ^2]  = v2  e”rt 


Therefore,  v must  satisfy 


2 -rt  .-Xt-  2 

v = / ( qv  + p)  e Ae  dt  = QV  +£ 

0 


Thus,  v is  one  of  the  two  positive  real  roots  Z and  Z'  given  by  (2.7) 
As  previously  shown,  Z < 1 and  Z 1 > 1.  Since  T0  > T it  follows  that 

□ 


v2  = E[e  rT2]  < E[e  *•]  = v.  Hence,  v = z; . 


,i  -rT'- 


Lemma  2.$  . Z = E[e  X],  i > 1. 


Proof  (induction  on  i) . Since  the  demand  process  is  memoryless, 


D 


= T'+T"  where  T'  and  T"  are  independent,  T'  = and  T"  = T.. 
Therefore. 


E[e'rTi+l]  = E[e  r(T,+T”)]  = E[e'rTl]  E[e‘rTi]  = ZZ1  = Zi+1 


□ 
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evaluated  at 


Thus,  Z1  is  simply  the  Laplace  transform  of  T^, 
the  discount  rate  r,  for  i > 0.  In  the  absence  of  intervening  permanent 
expansion,  the  excess  demand  (-k)  tracks  demand.  Hence,  if  there  is 
no  intervening  permanent  expansion,  T is  also  the  first  time  that  the 
spares  level  decreases  an  amount  i > 0 below  the  initial  spares  level. 
With  this  interpretation  in  mind,  the  alternative  probabilistic  proof  of 
Theorem  2.2  can  now  be  given. 


Proof  of  Theorem  2.2  (alternate).  Suppose  that  the  initial  spares  level 
is  k > 0.  Permanent  capacity  expansions  do  not  occur  while  the  spares 
level  is  positive.  Furthermore,  no  operating  costs  are  incurred  when 
the  spares  level  is  positive.  Therefore,  no  costs  can  accrue  until  time 
Tfc,  when  the  spares  level  first  becomes  nonpositive  (0).  Since  the  model 
is  time- stationary  and  the  demand  process  is  memory  less, 

ck(x,K)  = E[costs  during  [0,TR)  + costs  during  [T^co)] 

= E[  0 + e'rTk  C0(X,K) ] 

= E[e'rTk]  CQ(X,K) 

= Zk  CQ(X,K),  k > 0 . □ 


Let 


7 _ 1 - v/l  - 'iQP 

23 


(2.8) 
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The  following  lemma  summarizes  results  analogous  to  those  previously 
derived  for  Z. 

Lemma  2.6. 

( i)  Z is  real  and  a < Z < Q4-0, 

(ii)  ZL  = E[elT-i])  i > 0, 

( iii)  Z = p-IQZ. 

Proof.  (i)  follows  in  a manner  analogous  to  the  first  proof  of  Theorem  2.2 , 

~ 2 
as  Z is  the  smallest  root  of  the  quadratic  expression  |3v  - v + OC  = 0. 

( ii)  follows  in  a manner  analogous  to  the  proofs  of  Lemmas  2.4  and  2.5 

(iii)  follows  by  inspection  of  (2.6a)  and  (2.8).  □ 

2.4.  Recursive  Computation  of  the  Functionals  (Cq(-,K)) 

When  the  spares  level  is  nonpositive,  the  costs  behave  according 
to  (2.4)  and  (2.5).  Since  X > K,  it  follows  from  Theorem  2.2  that  (2.5) 
may  be  rewritten  as 

C_K(X,X)  = f_k(k)  + qc_k+1(x,k)  * p(tfe(K)  + g(x)  + ZX-K  C0(X,K)).  (2.9) 

Also,  C^(X,K)  = ZCq(X,K),  so  that  the  first  equation  of  (2.4)  can  be 
rewritten  as 


Cq(X,K)  = ( 1-QZ)  F0(K)  + e(l-QZ)-1  C_t(X,K)  . 
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Let  X = X and  K = K be  fixed.  Denote  C = C.  (X  K)  and 
Fk  = Fk(X),  0 > k > -K.  Then  (2.9)  and  (2.4)  yield  a square  (K+l) 
nonsingular  linear  system,  as  depicted  in  Figure  2.2.  Solving  this 

A A A 

system  yields  the  costs  (Ck)  for  the  fixed  values  X and  K.  By 
Corollary  2.3,  it  suffices  to  compute  only  CQ. 

Obviously,  it  is  impractical  to  solve  such  systems  for  all 

/\  /\ 

possible  values  (X,K),  so  some  sort  of  parameterization  is  necessary. 
Suppose  that  the  last  equation  in  the  system  of  Figure  2.2  is  replaced 
by  an  equation 


A /V  A A 

'^-K+l  + C-K  ’ gCd  = F 


+ £0e(K)  , 


where  is  a dummy  variable.  Then  the  resulting  linear  system  is 

/N  /N  /\ 

independent  of  the  value  X and  solving  for  CQ  in  terms  of  will 

A A 

/\  /N  /N/X  /\  X»K  ZV 

yield  CQ  = a + bC^.  Making  the  substitution  = g(X)  - Z CQ  then 

Z\ 

Z\  /\  /A  /\  ZN  X-K  -1 

yields  = (a  + bg(X))  ( 1 - bZ  ) . Since  this  is  true  for  any 


ZV  ZS  Z\ 


,X-Kx - 1 


value  X chosen^  it  follows  that  C^(X^K)  = (a  + bg(X))  (1-  bZ*  )"  y 

/\ 

for  all  X > K.  Hence,  using  the  dummy  substitution,  the  system  need 

A 

only  be  solved  once  for  a given  value  K to  obtain  the  functional 

ZN  zs 

Cg(*,K),  which  can  then  be  minimized  over  [K,oo)  . 

Referring  to  Figure  2.2,  notice  that  the  respective  linear  systems 

A A 

for  K = K and  K = K+l  would  be  very  similar.  In  fact,  the  coefficients 

A A 

in  the  first  K rows  and  K columns  would  be  identical.  Hence,  it  seems 
reasonable  to  expect  that  solving  for  the  functional  CQ( • ,k)  should 
aid  in  solving  for  the  functional  C^( • , K+l).  That  is,  it  should  be 
possible  to  recursively  solve  for  the  functionals  Cq(*,K)  for 


FIGURE  2.2:  Cost  Transition  Equations,  Model  I (X 


: 


K = 1,  2,  •••  • ,rhe  derivation  of  that  recursion  is  the  subject  of 

this  section. 


Definition  2,j.  Let  QI  and  0 satisfying  (2.2)  be  given.  Define  the 
transformation  : M+  x [0,1]  ->B+  by  W1(e,6)  = (1-06)"1  0e.  Define 
the  transformation  : [0,  1]  ->  [0,  1]  by  W2(5)  = a( 1 - 06) -1.  Define 
the  transformation  W : x [0,1]  ->R+  x [0,1]  by  W(e,6)  = (W^  e,6)  ,V*2(  6)  ) . 

The  following  lemma  verifies  that  the  above  transformations  are 
well-defined. 


Lemma  2.7.  ^^(0,6)  = 0,  W^e^)  >0  and  0 < W2(5)  < 1,  for  all 

€ > 0,  ° < 6 < 1. 

Proof.  By  Lemma  2. 1,  0 < 6 < 1 =>  1-06  >l-0>a>O=>O<(  1-06)  "1  < a-1 
0 < a(  1-06)  < 1 =>  0 < Wg( 6)  < 1.  Note  that  W^(e>6)  = ( 1-06) ”1  0e 

= a X0  W2(6)e.  Therefore,  0 < 6 < 1 =»  W1(0,6)  = 0 and  W1(e,6)  > 0, 
for  e > 0.  □ 

Since  W is  an  into  mapping,  the  composition  mapping  w”*  is 
well  defined:  W°(e,5)  = (e,6)  and  Wm(e,6)  = W o Wm"1(e,6),  m > 1. 
Similarly,  the  composition  mapping  W™  is  well-defined:  W°(6)  = 6 

and  W2(6)  = W2^_1(&),  m > 1.  By  the  previous  lemma,  the  following 
transformations  are  also  well-defined. 


A 
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Definition  2.4.  Let  a and  3 satisfying  (2.2)  be  given.  Let 

=IR  X [0, 1]  x (0,  1,  . . . , M]  xlR^.  For  M = 1,  2,  . ..,  define  the 
transformations  S : 5^  -» 


by 


S(p,5,i,V) 


(P,  6,  V)  , 

(a-1  w2(B)  (3P+V.),  w2(6),  i-i,  v), 


i = 0 
1 < i < M 


where  V = (Vp  Vg,  . ..,  V ) £ 1RM. 


The  above  definition  defines  a separate  transformation  for  each 
positive  integer  M.  To  avoid  possible  ambiguities,  it  will  be  assumed 
that  dim  V = M implies  which  transformation  is  being  used.  The  composi- 
tion mappings  Sm  are  well-defined  for  all  m > 0:  S°(p,5, i,V)  = (P,&,i,V) 

and  Sm(p,5,i,V)  = Sm  ^ o S(p,6,i,V).  S™( p,5, i,V)  will  denote  the  jth 

component  of  Sm( P,&, i,V) , j = 1,  2,  3. 


Lemma  2.8.  For  all  p £1R,  & £ [0, 1]  and  vectors  V, 

sJPjM*7)  = w“(&)  > 0 < m < i < dim  V , 

S®(p,&,i,V)  = i-m  , 0 < m < i < dim  V . 

Proof  (induction).  Let  P £IH,  & £ [0, 1],  M £ {1,  2,  ...)  and 
M 

V £1R  be  arbitrarily  chosen.  For  m = 0,  the  lemma  is  trivial  and 
for  m = 1,  the  lemma  follows  from  Definition  2.4  since  i > m = 1 and 

i < dim  V = M.  Assume,  in  general,  that  the  lemma  is  true  for  some 
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m- 1 < M.  Denote  P'  - (P,&,i,V).  Let  i > m and  i < M.  By  the 

induction  hypothesis. 


s,V,M,v)  = s « 


S(P',  V^-L(5),  i-(m-l),  V)  . 


Since  i > m ^ i-(m-l)  > 1,  it  follows  from  Definition  2.4  that 

S”(p,M,V)  = w2  ° S^P.M,^  = W2  ° w2_1(5)  = W?&)  and  s^(  P, &,  i,V) 
= (P>^,i,V)  - 1 = i-(m-l)-l  = i-m.  Hence,  the  lemma  follows  by 

induction  on  m.  n 


The  transformations  S have  two  important  properties  that  permit 
a simple  recursive  solution  of  the  problem  at  hand.  These  properties 
are  given  in  the  next  two  lemmas. 

Lemma  2.9.  Let  V = (Vp  ...,  VM)  elRM  and  let  V'  = (0,V)  €1RM+1.  Then 


SlK6,i,V)  = S®(p,&,i-l,V), 


l<m<i<M+l 


Proof  (induction).  For  m = 1,  let  i > m.  Then  i > 1,  so  ^ 

and  i-1  > 1.  Hence,  by  Definition  2.4, 

slKM,V')  = a'1  »2(6)  Op+v^) 

= a'1  w2(6)  (pp+v^)  = S^( P,6, i-  1,V)  . 


iAsr. 


— - — I 1 1 


JVf 


Pro°^  (induction).  For  m = 0,  B(0,b,M)  = (unique),  where  e.  denotes 

M ^ 

the  jth  unit  vector  of  IR  . In  general,  assume  the  lemma  is  true  for 

some  m-1  > 0,  m-l<M-l.  Then  S^' L( VM,b,M- 1,V)  = B(m- l,b,M) • V. 

Denote  b'  = s™"  \ v) • Since  m-1  <M-1,  Lemma  2.8  gives 

s™  ^V)  = M-m  > 1-  Thus,  by  Definition  2.!+, 

SI(V6>M_1>V)  = a_1  W2(6')  + VM_J 

= B(m,b,M)-V  , 

where  (uniquely) 

B(m,b,M)  = a'1  W (&')  (£}B(m- 1,5  M)  + eM  ) . 

Thus^  the  lemma  follows  by  induction  on  m.  □ 

Denote  F(K)  = (Fq(K),  F_1(K),  . ..,  F_K(K)).  Then  F(K)  £IRK+1 
and  Fk(K)  is  the  -(k-l)th  component  of  F(K) , 0 > k > -K. 

Theorem  2.11. 

ck(X,K)  = Pk(K)  + Br(K)  Ck+1(X,K)  + ek(K)  (g(X)  + p&( K)  + ZX-KCq(X,K))  , 

x > K > °,  0 > k > -K  (2.10) 

where 
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and 


Uk(K),  &k(K))  = ^(3,00  , 


(2.10a) 


Pk(K)  = Sx+  (F_r(K),  a,  K,  F(K) ) . (2.10b) 

Froof  (induction).  Denote  U(X,K)  = g(X)  + 0e(K)  + ZX_K  CQ(X,K). 

For  k = -K,  W°(pfa)  = (p,a)  and  sJ(F_R(K),  a,  K,  F(K) ) = F_r(K), 
so  (2.10)  agrees  with  (2.9).  Hence^  the  theorem  is  true  for  k = -K. 

In  general^  suppose  that  the  theorem  is  true  for  some  k-1,  -K  < k-1  < -1. 
Then, 

WM)  = pk.i<K)  + \.i(K)  ck(x>K>  + ek-i(K)  U(X>K>  * (2. 11) 


By  (2.V), 


ck(x,K)  = Fk(K)  + 0Ck+1(X,K)  + fiC^^K)  . 


(2.12) 


Substituting  (2.11)  into  (2.12)  gives 


Ck(X’K>  = Pk(K)  + VK>  Ck+1(X’K)  + ek(K)  U<X>K>  > 


where 


€k(K)  = (l-P&^ifK))'1  fJe^K)  = W^e^K),  ^00)  , 
&k(K)  = 0£(  1-P&k_ x( K) ) ” 1 = Wg(&k_1(K))  , 


Pk(K)  = (l-e&j^K))"1  OPk-1(K)  + Fk(K)) 
= a~l  W2(Sk-l^K))  ^pk.!(K)  + F (K)) 


Hence,  using  the  induction  hypothesis, 

<VK><  Vc»  * “(Vi<K>>  Vi<*» 

- “ • "K’hk‘  l(  B,«)  . wK*k(Bja)  > 

and 

(pk(R),  6k(R),  -k,  F(K) ) = S(Pk_l(K),  &k  l(K),  -(Id),  F(j 

= S o SK+k-L(F_K(K),  Of,  K,  F(K)) 

= sK+V_k(K),  a,  K,  F(K) ) . 


The  positivity  of  (ejK),  6JK))  follow  from  Lemma  2.7.  Thus,  the 
theorem  follows  by  induction  on  k.  Q 

Theorem  2.11  provides  an  immediate  recursion  for  the  coefficients 
e0(0  and  &Q(.). 

Corollary  2. 12. 

(i)  eQ(0)  = 0;  e0(K+l)  = (1  - g&^K))'1  £k0(K),  K > 0, 

(ii)  &Q(0)  = a;  &0(K+i)  = (1  - B6q( K) ) ~ 1 a,  K > 0. 
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Ir.2£.f-  (e0(0),  60(0))  = WU(0,a)  = (0,a).  For  K > 0, 


(e0(K+1)>  &o(K+1))  = = W o W^a)  = W(e0(K),  &0(K)) 

= ((i-eSoCK))"1  6e0(K),  (1  - P&^K))'1  a)  . 


Lemma  2. 13. 


1 - &0(K)Z 

eo(K> 


K > 0 . 


Proof  (induction).  Recall  that  Z satisfies  -OZ2+Z-0  = 0,  so 
1-QZ  = 0Z  \ Using  Corollary  2.  12f 


~ V°)Z  z»  l-QZ  BZ'1 

0(0)  * B ■ — ■ Z 


Suppose^  in  general  that  the  lemma  is  true  for  K > 0.  Then  using 
- 1 2 

0 (Z-QZ  ) = 1 and  Corollary  2.  12^ 


1 - &q(K)Z 

6oW 


1 - &0(K>Z  ZK+1 

Zeo(R) 


0~  L( Z-OZ2)  - 60(K)Z  ^K+1  1 - QZ  - 06o(K)  K+1 

” Ze0(K)  ' Pe0(K)  2 

1 - (i  - es^K))-1  OZ  zK+1  1 - &0(K+1)Z  zK+1 
(1-06JK))-1  0e  ( K)  = e0(K+1)  Z 


r 


H 


i i 


Using  the  proceding  lemma^  the  functional  form  of  Cq(*)K)  can 
not.  be  obtained. 


Theorem  2. 14. 


c (XK)  - OT  t--g(X).  zK 
c0(x,k;  - -1  x » 


X > K > 0 


z - z 


(2.11) 


where 


po(K> 

*<K>  =Pe(R)  +7JkT 


(2.11a) 


Proof.  By  Theorem  2.11^ 


X-K. 


Cq(X,K)  = Pq(K)  + 6q(K)  C][(X^ K)  + €o(K)(0e(K)+g(X)4€o(X,K)ZA-*) 


By  Theorem  2.2 f C^(X^K)  = ZCq(X>K).  Substituting  this  result  into  the 
above  expression  and  collecting  terms  yields 


yx,K)  = 


(0e(K)  + P0(K)/e0(K))  + g(X) 


(1-60(K)Z)  ZK/6q(K)  - ZX 


whichj  by  Lemma  2. is  (2.11) 


□ 


Recall  that  Corollary  2.11  gives  a simple  recursion  for  £q(*). 
Hence?  by  (2.113)^  a similar  recursion  for  Pq( ' ) will  provide  the 
means  for  recursively  computing  the  functionals  C^( • ,K) , K.  = 0^  1^  2f 


The  next  theorem  provides  the  necessary  recursion  for  p^( • ) 


1 
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Theorem  2.  H. 


Pq(K)  = B(K)-P(K)  , K>0, 


(2.12) 


where 


b(o)  = (i)  em1 


B(K+1)  = a_1  P50(K+1)  (fTl,  B(K) ) €BK+2  , K > 0 . 


(2.12a) 


Proof.  By  Lemma  2.10  and  (2.10b),  the  vectors  B(K)  = B(K,a,K+l)  £IRK+F 
exist  and  are  unique,  K > 0.  By  (2.10b),  pQ(0)  = FQ(0),  so  B(0)  = (1). 
Recall  that  F(K)£lRK+1  and  F(K+1)  £ ]RK+2.  Denote 

F(K+1)  = (f_1(k+i)  , ...,  F_  (K*l))  £ mK+1  . 


By  (2.10b)  and  Lemma  2.9, 


P-i(K+l)  = Vf-(k+i)(k+1)>  a;  K+1,  F(R+1)) 
= si(f-(k+i)(k+1>>  a>  K>  F(K+1))  • 


Since  F(K+1)  6 1R  * } Lemma  2.  10  gives 


P_X(K+1)  = B(K,a,K+l)  • F(K+1)  = B(K)  • F(K+1)  . 


By  Lemma  2.8  and  (2.10a), 


S2^F.(K+1)^K+1^  a>  K+1>  F(R+1))  = w|(a)  = Sq(K)  . 

Hence, 

sK(f-(k+i)<k+1)>  a>  K+1,  p(K+i)) 

* (B(K)  • F(K+1),  5q(K),  1,  F(K+1) ) . 

Thus,  by  (2.10a)  and  Corollary  2.12, 

P0(K+1)  = S1(B(K)-F(K+1),  Bq(K),  1,  F(K+1) ) 

= a_1  W2(60(K))  OB(K).f(K+l)  + Fq(  K+l) ) 

= ex'  1 e&0(K+l)  ( B(K)  •F(K+1)  + fT1  Fq(K+1)) 

= a'1  f»0(K+l)  O'1,  B(K) ) • (F0(K+1),  f,K+l)) 

= Q1’1  3&0(K+1)  (0_1,  B(K) ) • F(K+1)  . 

Denote  0(K)  = of1  e&Q(K);  then,  BQ(K)  = p"1  O0(K).  The  follow- 
^8  theorem  summ3ri.zes  the  required  recursions 

Theorem  2.15. 

(i)  6(0)  = 0;  e(K+l)  = 0(1  - ae(K))'1,  K > 0, 

(ii)  eQ(0)  = Q;  e0(K+l)  = 0(K+1)  eQ(K) , K > 0, 

(ill)  B(  0)  = (1);  B (K+l)  = 0(K+1)  (0_1,  B(K) ) , K > 0. 

6? 


Ill  I III  I ■ I .1 


_Proof.  Note  that  0(1-  0bQ(K))  1 = 0(  1 - 00_  la6<  K) ) " 1 = 0(1-  08(K))'1. 
By  Corollary  2. 12,  0(0)  = cfl0&  (0)  = 0 and 

0(K+l)  , a_10&o(K+i)  = a'W(i  - e&0(K))_1 

= 0(1  - 06o(K))_1  = 0(1  - ae(K))*1,  for  K > 1 ; 

this  proves  ( i) . By  ( i)  and  Corollary  2. 12,  eQ(0)  = 0 and 

e0(K+l)  = (1  - 0&q(K))-1  0eo(K)  = (l-a0(K))_1  eQ( K)  = e(K+l)  e0(K)  ; 

this  proves  ( ii) . By  Theorem  2. 15 

B(K+1)  = a"106()(K+ 1)  (0_1,  B( K) ) = 0(K+1)  (0*1,  B( K) ) , 

which  proves  ( iii) . □ 

By  Theorem  2. 14,  the  functional  C^( • ,K)  is  completely  determined 
by  the  coefficient  0(K).  Hence,  a recursive  procedure  for  computing  j6(  • 
is  a recursive  procedure  for  computing  the  functionals.  Using  Theorems 
2.14,  2.15  and  2. 16,  the  desired  algorithm  can  now  be  stated. 
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Algorithm  Al  (Compute  j6(K),  K = 0,  1,  2,  . . . , K) . 

(1)  K «-  0,  0 «-  p,  e <-  0,  B *-  ( 1)  . 

(2)  p *-B-F(K),  0(K)  <-  e” *P  + 0JK) . 

(3)  K <—  K+l,  e <-  0(  l-3e) ' 1)  e <-  6ef  B 

(*0  If  K > K,  stop;  otherwise,  go  to  step  (2). 

Although  the  recursion  for  PQ( • ) can  be  stated  in  a number  of 
alternative  ways,  the  vector  dot  product  is,  in  general,  unavoidable. 

However,  when  the  expected  incremental  costs  are  both  proportional  and 

constant  in  K,  then  the  dot  product  can  be  eliminated,  as  will  now  be 
shown. 

Denote  AFk(K+l)  = F^K+l)  - FjK),  0 > k > -K,  K > 0.  Denote 
AF(K+1)  = (AF0(K+1),  ...,  AF_k(K+1))  £KK+1.  Then , F^lC+l)  = FjK) 

+ AFk(K+l),  k = 0,  -1,  ...,  -K.  Hence,  F(K+1)  = (Fq(K+1),  F(K) ) +(0,AF(K+1) ) , 

K > 0. 


Lemma  2. 17 


P0(K+1)  = 0(K+1)  O'1  Fq(K+1)  + B(K).AF(K+1)  + PQ(K))  , K > 0.  (2.13) 


Proof.  By  Theorem  2. 1 6 


p0(k+1)  = B(K+1).F(K+1)  = &(K)  O-1,B(K)).((F0(K+l),F(K))  + (0  ,AF(K+1)) 

= 0(K)  (0'1  F0(K+1)  + B(K).F(K)  + B(K)  *AF(K+1) ) 

= 0(R)  O"1  Fq(K+1)  + P0(K)  + B(  K)  *AF( K+l) ) . □ 


Recall  that,  when  the  incremental  costs  are  constant  in  K,  it 
can  be  assumed  that  Fk(K)  = F 0 > k > -K,  K > 0 (by  the  introduction 
of  the  function  0e(K)).  If  the  incremental  costs  are  also  proportional 
(see  Example  2.3),  then  AF(K)  = AF  1^  where  AF  is  a scalar  constant 
and  1^  denotes  the  K-vector  of  all  l's.  In  this  case,  (2.13)  reduces 
to 

P0(K+1)  = 0(K+1)  (0_1FO  + AF  b(K)  + P0(K)),  K > 0 , (2.14 

where  b(0)  = 1,  and 

b(K)  = B(K)-1K+1  = e(K)  (ff1,  B(K-1)).1K+1 

= e(K)  (0_1  + B(K-1).1K)  = 0(K)  (3'1  + b(K-l)),  K > 1 . (2.14 

ihus,  the  vector  dot  product  is  eliminated  when  the  expected  incremental 
costs  are  both  proportional  and  constant  in  K,  giving  the  following 
algorithm. 

Algorithm  Ag  (Compute  0(K),  K = 0,  1,  2,  . . . , K,  when  F's  are 
Proportional  and  Constant  in  K) . 

( 1)  K <-  0,  0 B,  P «-  Fq,  6 «-  p,  b «-  1. 

(2)  j6(K)  <-e_1p  + 0e(K). 

(3)  K e-  K+  1,  e <-B(  1-pe)"1,  £<-&€,  P «-  0(  P”  1Fq  + AFb  + P)  , b 4-  0(  &"  X+b) 

(4)  If  K > K,  stop;  otherwise,  go  t?  step  (2). 


Given  the  recursions  of  this  section,  it  is  possible  to  derive 
closed-form  expressions  for  the  coefficients  0(*)  defining  the 
functionals  C^(*,K),  K _>  0.  The  derivation  is  accomplished  by  studying 
the  determinants  of  some  very  special  tridiagonal  submatrices  in  the 
linear  system  depicted  by  Figure  2.2  [16].  However,  from  a computation 
aspect,  these  expressions  are  of  limited  value,  since  they  do  not 
apparently  provide  any  efficiencies  over  the  recursive  algorithms  given 
here. 


2.5.  Assumption  2.5  with  Economic  Interpretation 


This  section  discusses  the  final  assumption  of  Model  I: 


Assumption  2.3.  The  coefficients  0(-)  of  (2.11)  are  nonnegative  and 
nondecreasing. 

The  form  of  the  functional  Cq(*,K)  can  be  anticipated  through 
a probabilistic  argument  that  provides  an  economic  interpretation  of  the 
coefficient  0(K).  Given  the  economic  interpretation,  it  is  demonstrated 
that  the  above  assumption  in  essence  states  that  the  costs  of  temporary 
facilities  exceed  the  permanent  facility  operating  charges  that  would 
otherwise  be  incurred. 

Let  G^K)  denote  the  expected  discounted  costs  until  the  first 
expansion,  starting  from  spares  level  k > 0,  when  K is  the  limit  on 
temporary  facility  usage. 


Proof.  Starting  from  spares  level  k > 0,  no  costs  are  incurred  until 
the  spares  level  first  becomes  nonpositive;  the  time  of  this  event  is 
given  by  Tk  end  EU'^)  . zk  by  Le™a  2.i».  The  expected  costs  beyond 
time  Tk,  until  the  first  expansion,  are  given  by  Gq(K)  (discounted 
relative  to  Tk) . Hence 


Gk(K)  - E[Gq(K)  e k]  = G0(K)  E[e  “k]  = GQ(K)  Z 


-rT,. 


When  an  expansion  occurs , the  spares  level  changes  to  X-K  > 0.  Hence, 
the  inter- expansion  costs  are  given  by  G^  ^(k)  = ZX_X  G (K) 


Theorem  2.19. 


C0(X,K)  = 


Gq(K)  Z-(R+1)  + g(x) 


Z'1  - zx 


X > K > 0 


(2.16) 


Proof;  Starting  from  spares  level  zero,  the  time  until  the  first  expansion 
is  given  by  At  time  TR+1,  an  expansion  of  size  X+l  occurs  at 

cost  g(X)  and  the  expected  costs  thereafter  are  C (X  K)  = ZX_K  C (X  K) 
(discounted  relative  to  ) . Hence, 


C0(X,K)  = Gq(K)  + E[(g(X)  + ZX’K  CQ(X,K))  e'  K+1] 
= Gq(K)  + (g(X)  + ZX'K  Cq(X,K))  E[e'rTK+1] 
= Gq(K)  + (g(X)  + ZX'K  CQ(X,K))  ZK+1 
= Gq(K)  + g(X)  ZK+1  + ZX+1  Cq(X,K) 

G (K)  + g( X)  ZK+1  G (K)  Z^K+1)  + g(X) 

n / v u 


C0(X,K)  = 


i - r 


z-1  - zx 


Corollary  2 . 20 . 


0(K)  = Gq(K)  Z^K+1)  , 


K > 0 


(2.17) 


Proof.  This  follows  by  comparing  (2.16)  and  (2.11),  each  of  which  is 
satisfied  for  arbitrary  g and  X > K > 0 . □ 


Rearranging  (2.17)  gives  Gq(K)  = 0(K)  ZK+1  = 0(K)  E[e’rTK+1]. 
Similarly,  using  (2.15)  gives  G^(K)  = 0(K)  zK+k+1  = 0(K)  E[e'r^+k+1], 
k > 0.  When  the  initial  spares  level  is  k > 0,  the  time  of  the  first 
expansion  is  Thus,  in  an  expectation  sense,  0(K)  is  the 

equivalent  lump  sum  payment  that  would  be  incurred  at  the  first  permanent 
expansion  time  in  lieu  of  the  incremental  charges  over  the  time  span 
previous  to  the  expansion.  Since  the  expected  inter- expansion  costs  are 
Gjj_^(K)  = 0(K)  E[e  X+lj  and  tfte  inter-expansion  time  intervals  are 
equal  (in  distribution)  to  a similar  lump  sum  interpretation  of 

0(K)  can  be  given  with  regard  to  inter- expansion  costs. 
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Corollary  2.21.  For  K > 0 and  k > 0, 

(i)  0(K)  > 0^Gk(K)  > 0, 

(ii)  0(K+1)  >0(K)  ^Gr(K+1)  > ZGr(K), 

(iii)  Gk(K+l)  > Gk(K)  =>P(K+1)  >0(K). 


Proof.  ( i)  follows  directly  from  (2.17)  and  (2.15).  ( ii)  follows  from 

(2.17)  as  p(K+l)  - 0(K)  = (GQ(K+1)  - Gq(K)Z)  z'(K+1)  = (Gfc(K+l) 

- Gk(K)Z)  Z \ ( iii)  follows  from  ( ii)  since  0 < Z < 1.  □ 


By  ( i)  and  (iii),  the  coefficients  $(•)  will  be  nonnegative  and 
nondecreasing  if  the  expected  inter-expansion  costs  are  likewise.  Since 
the  only  non-expansion  costs  are  the  expected  incremental  costs  incurred 
when  temporary  facilities  are  in  use,  ( i)  states  that  p(K)  will  be 
nonnegative  if,  in  aggregate  expectation,  the  costs  of  the  temporary 
facilities  used  exceed  the  permanent  facility  operating  charges  that 
would  otherwise  be  incurred.  With  regard  to  (iii),  assume  that  the 
initial  spares  level  is  kQ  > 0.  Using  temporary  facilities  limit  K, 
the  time  of  the  first  expansion  is  TK+ko+]/  Using  temporary  facilities 
limit  K+l,  the  time  of  the  first  expansion  is  T^+k  +^.  Hence,  Gk(K) 
represents  the  expected  costs  over  the  time  interval  [0,  and 

Gk(K+l)  represents  the  expected  costs  over  the  time  interval 
[0,  TK+ko+2^  ~ ^K+k  +1^  (with  probability  1).  It  is  difficult  to 
imagine  practical  situations  where  the  expected  costs  over  the  time 
interval  [0,  T^+k  + using  limit  K could  exceed  those  over  the  greater 
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time  interval  [0,  T^^  usi-n8  limit  K+l. 

indicates  that  there  is  really  no  optimization 


In  fact,  the  next  lemma 
problem  if  0(K+1)  <0(K). 


Lemma  2.22.  Let  K > 0.  If  0(K+1)  <0(K),  then  CQ(X,  K+l)  <CQ(X,K), 
X > K+l. 


Proof.  By  (2. 11),  for  X > K+l, 


c„ (x,k+d  = ± jw  < fty  sw 

u z'1  - tr  z-i  - zx 


c0(x,k)  . 


□ 


Thus,  if  K and  K+l  are  each  potential  limits  on  temporary 
facility  usage,  it  suffices  to  compare  CQ(K,K)  with  {CQ(X,K+1), 

X > K+lj.  Thus,  from  an  optimization  standpoint,  Assumption  2.3  simply 
eliminates  trivialities. 


2 .2 * * *  6 • Summary  and  Statement  of  the  Expansion  Size  Optimization  Problem 
Let  k denote  the  set  of  feasible  values  for  K.  To  illustrate, 

< = {0,  1,  2,  ...,  K]  would  be  typical  of  Examples  2. 1 - 2.4,  where  K 

denotes  an  upper  bound  on  temporary  facility  usage  based  on  physical 

considerations.  Similarly,  k = {0,  L,  2L,  ...,  SL)  for  the  restricted 
modular  Example  2.5,  where  N denotes  an  upper  bound  on  temporary 
module  usage  based  on  physical  constraints. 

Given  K £ k and  an  initial  spares  level  > 0,  an  optimal 
permanent  expansion  size  X*(K)  + 1 for  Model  I is  given  by 

75 


i 


f 


l 


CQ(X*(K),  K)  = min(CQ(X,K)  : integer  X > K)  , 


(2.18) 


where 


CQ(X,K)  = i '+  zK  , 


X > K } 


(2.19) 


&(')  is  nonnegative  and  nondecreasing  over  k : 


(2.20) 


0 < Z = t/1-4c*£_  ^ ^ given  a and  (3  satisfying 


(2.2) 


( 02,  - Z + £ = 0); 


(2.21) 


Ck(X,K)  = Zk  C0(X,K),  k>0;  (2.22) 

and 

- rT 

Zk  = E[e  k]  , k > 1 . (2.23) 

(2.19)  follows  from  Theorem  2. 14,  (2.20)  is  Assumption  2.3  (discussed  at 
length  in  the  previous  section),  (2.21)  and  (2.22)  follows  from  Theorem 
2.2,  and  (2.23)  follows  from  Lemma  2.5.  The  fact  that  X*(K)  is  given 
by  (2.18),  independent  of  k^  > 0,  follows  from  Corollary  2.3. 

Algorithms  A^  and  Ag,  presented  in  Section  2.4,  provide  simple 
recursive  procedures  for  determining  the  coefficients  0(*)  over  k. 

Once  these  coefficients  are  known,  the  task  of  determining  X*(*)  over  k 
becomes  the  series  of  single-variable  minimization  problems  given  by 
(2.18)  under  conditions  (2.19)  - (2.23).  This  optimization  problem  is 
treated  in  [17].  By  Corollary  2.3,  the  optimal  temporary  facilities 
usage  limit  K*,  and  the  associated  optimal  expansion  size  X*+l,  are 
then  given  by 


CQ(X*,  K*)  = min{C0(X*(K),  K)  : K € «]  . 
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facilities  represent  the  extraordinary  measures  taken  in  order  to  serve 
excess  demand  (prior  to  an  expansion  of  permenent  facilities).  Examples 
of  temporary  facilities  include  backlogging,  overloading  and  jobletting.  This 
paper  considers  non-modular  temporary  facilities,  the  costs  of  which  depend 
only  on  the  amount  of  excess  demand  currently  being  served.  A companion 
paper  [17]  treats  the  case  of  modular  temporary  facilities. 

For  a given  limit  (K)  on  temporary  facility  usage,  the  form  of  the 
expected  discounted  cost  functional,  parameterized  in  the  expansion  size 
(X+l),  is  derived.  Recursions  are  given  for  determining  these  functionals 
over  all  feasible  values  for  K.N  A sequel  paper  [18]  treats  the  problem  of 
minimizing  the  functionals  in  order  to  find  optimal  expansion  sizes. 
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